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ABSTRACT 

Helioseismology provides a powerful tool to explore the deep interior of the Sun. Measure- 
ments of solar interior quantities are provided with unprecedented accuracy: for example, the 
adiabatic sound speed c can be inferred with an accuracy of a few parts in 10 . This has be- 
come a serious challenge to theoretical models of the Sun. Therefore, we have undertaken a 
self-consistent, systematic study of sources of uncertainties in the standard solar model, which 
must be understood before the helioseismic observations can be used as constraints on the theory. 
This paper focusses on our own current calculations, but is also a review paper summarizing the 
latest calculations of other authors. We find that the largest uncertainty in the sound speed c in 
the solar interior, namely, 3 parts in 10 3 , arises from uncertainties in the observed photospheric 
abundances of the elements: C, N, O, and Ne have uncertainties of ~ 15%, leading to an uncer- 
tainty of ~ 10% in the photospheric Z/X ratio. Uncertainties of 1 part in 10 3 in the sound speed c 
arise, in each case, from (1) the ~ 4% uncertainty in the OPAL opacities, (2) the ~ 5% uncer- 
tainty in the basic pp nuclear reaction rate, (3) the ~ 15% uncertainty in the diffusion constants 
for the gravitational settling of helium, and (4) the ~ 50% uncertainties in diffusion constants 
for the heavier elements. (Other investigators have shown that similar uncertainties arise from 
uncertainties in the interior equation of state and in rotation-induced turbulent mixing.) In the 
convective envelope only, uncertainties in c of order 1 part in 10 3 arise from the uncertainty of a 
few parts in 10 4 in the solar radius, and from uncertainties in the low-temperature equation of 
state. Other current uncertainties, namely, in the solar age and luminosity, in nuclear rates other 
than the pp reaction, and in the low-temperature molecular opacities, have no significant effect on 
the quantities that can be inferred from helioseismic observations. Significant uncertainty in the 
convective envelope position R ce (of up to 3 times the observational uncertainty of ±0.001 Rq) 
arises only from uncertainties in Z/X, opacities, the pp rate, and helium diffusion constants; 
the envelope helium abundance Y e is significantly affected (±0.005) only by extreme variations 
in Z/X, opacities, or diffusion constants, and is always consistent with the "observed" range of 
helioseismically inferred Y e values. Our predicted pre-main-sequence solar lithium depletion is a 
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factor of ~ 20 (an order of magnitude larger than that predicted by earlier models that neglected 
gravitational settling and used older opacities), and is uncertain by a factor of 2. The predicted 
neutrino capture rate is uncertain by ~ 30% for the 37 C1 experiment and by ~ 3% for the 71 Ga 
experiments (not including uncertainties in the capture cross sections), while the 8 B neutrino 
flux is uncertain by ~ 30%. 

Subject headings: diffusion — neutrinos — Sun: abundances — Sun: helioseismology — Sun: 
interior 



1. Introduction 

Helioseismology provides a powerful tool to explore the deep interior of the Sun. Measurements of solar 
interior quantities are provided with unprecedented accuracy. The Michelson Doppler Imager (MDI) on the 
Solar and Hcliosphcric Observatory (SOHO) spacecraft (Rhodes et al. 1997), the Global Oscillation Network 
Group (GONG: see, e.g., Gough et al. 1996; Harvey et al. 1996), the Birmingham Solar Oscillation Network 
(BiSON: Chaplin et al. 1996), and the Low-Hnstrument (LOWL: Tomczyk et al. 1995) provide helioscismic 
frequency measurements with accuracies of a few parts in 10 5 . From these, the sound speed c throughout 
most of the solar interior can be inferred with an accuracy of a few parts in 10 4 , as can the adiabatic index T\; 
the density can be inferred with an accuracy of a few parts in 10 3 (Basu, Pinsonneault, & Bahcall 2000; 
Bahcall, Pinsonneault, & Basu 2001). The depth of the solar envelope convection can be measured with 
an accuracy of nearly a part in 10 3 (Basu & Antia 1997). The abundance of helium in the solar envelope 
can also be inferred, but this value depends on the solar equation of state (see, e.g., Basu & Antia 1995; 
Richard et al. 1998; Di Mauro et al. 2002); as discussed in § 3.3, this inferred helium abundance appears to be 
uncertain at the level of a few percent (i.e., an uncertainty of at least 0.005 in the helium mass fraction Y e ). 

The ultimate goal of our work was to explore systematically a wide range of solar models with relatively 
modest amounts of mass loss on the early main sequence, and to test their viability via helioseismological 
measurements (as well other observational constraints). This mass loss investigation is presented in our 
companion paper "Our Sun V" (Sackmann & Boothroyd 2002). Since the consequences of moderate early 
solar mass loss are expected to be small, it is essential to understand the consequences of uncertainties in 
the input physics and in the input parameters of the solar model. The physics inputs include the equation 
of state, opacities, nuclear rates, diffusion constants, the treatment of convection (including the possibility 
of overshoot), and the effects of rotation and mass loss. Observed solar parameters include the solar age and 
the present solar radius, luminosity, surface composition, and solar wind. A considerable amount of work 
has been published investigating many of the above effects, as will be discussed in § 3. However, some basic 
uncertainties still warrant further attention. Before proceeding to our mass loss work, we found it necessary 
to try to extend the investigations of the above uncertainties. In particular, the major consequences arising 
from uncertainties in the present observed solar surface Z/X ratio, in the basic p-p chain rate, and in the 
diffusion constants for gravitational settling of helium and the heavy elements, have not been addressed 
sufficiently in the recent (most accurate) work of other investigators. To determine the consequences of 
these uncertainties, we computed standard solar models with various values of the Z/X ratio, the p-p chain 
rate, and the diffusion constants (lying in the permitted ranges). To obtain a set of self-consistent results, 
we also investigated the consequences of a number of other uncertainties, namely, in other nuclear rates, in 
envelope opacities, in the equation of state, from different methods of handling interior opacities, and in the 
solar age, luminosity, and radius. It is the aim of this paper to present the effects of the above uncertainties 
on the run of the sound speed and density in the solar interior, and on the radius R ce of solar convection 
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and the solar envelope helium abundance Y e . We also present the effects on the pre-main-sequence solar 
lithium depiction (excluding rotation effects) and on the production of solar neutrinos. (Note that lithium 
depletion in a non-rotating standard solar model occurs entirely on the pre-main-sequence; this predicted 
depletion is significantly smaller than the observed solar lithium depletion factor. Rotational mixing on the 
main sequence is generally invoked to account for the remaining lithium depletion). 

2. Methods 

We computed a reference standard solar model using up-to-date physics plus the observed values for 
the solar parameters. Several dozen variant solar models were also computed, in which one of these "inputs" 
was varied within the allowed uncertainties. (In a few cases, more than one of the "inputs" was varied, or 
the size of the variation exceeded the size of the quoted uncertainty in order to get a better estimate of the 
sensitivity.) By comparing with the reference standard solar model, the sensitivity to the uncertainties in the 
"inputs" could be determined. The stellar evolution code used to compute these solar models is descended 
from that used earlier in our "Sun III" paper (Sackmann, Boothroyd, & Kraemer 1993; see also Boothroyd 
& Sackmann 1999), but has been extensively updated for improved accuracy, including provision for much 
finer zoning as well as up-to-date input physics. 

Equation of State: The reference standard solar model used the OPAL equation of state (Rogers, Swen- 
son, & Iglesias 1996) in the interior, and the MHD equation of state (Dappen et al. 1988) where it was 
designed to be accurate, namely, in the outer envelope at log p < — 2 (this corresponds to r > 0.94 Rq and 
log T < 5.5 in the present Sun). A version of the MHD equation of state computation program was kindly 
provided to us (D. Mihalas 1999, private communication), and minor modifications allowed computation of 
the MHD equation of state for various hydrogen abundances and metallicities down to pre-main-sequence 
photosphcric temperatures and even below (A. I. Boothroyd, in preparation). In both cases, the equation of 
state was interpolated in metallicity Z as well as in hydrogen abundance X, temperature T, and density p, 
in order to take into account metallicity variations due to diffusion and nuclear burning. Tables for the 
MHD equation of state were computed on the same {X, Z, T, p) grid as the OPAL tables, with the following 
exceptions: for MHD, the upper p limit was slightly higher than for OPAL at low temperatures, but lower at 
high temperatures (since it never exceeded log p = —1.5); below logT — 3.7, intervals of A log T — 0.01 were 
used for the MHD T-grid, the same as the average log T interval in the OPAL grid for 0.005 < Tq < 0.006 
(3.7 < log T < 3.78); and for 0.006 < T 6 < 0.05 (3.78 < log T < 4.7), extra T-grid points were added to 
double the density of the grid as a function of temperature in this region. By computing some MHD tables 
at intermediate grid-values, the interpolation accuracy in each of {X, Z, T, p} was tested. In the outer solar 
convective zone (r > 0.95 Rq), the fractional error in the pressure P was < 2 x 10~ 4 , and for derivative 
quantities such as C v , \t, Xpi an d Pi, the fractional error was < 2 x 10~ 3 ; deeper in the Sun (r < 0.95 Rq), 
the interpolation errors were nearly an order of magnitude smaller. The interpolation accuracy is thus several 
times better than the quoted precision of the OPAL equation of state in the corresponding regions (Rogers 
et al. 1996). 

Variant cases were tested where the OPAL equation of state was used down to log T = 4 or even all the 
way down to its lower limit of validity at log T = 3.7 — in the latter case, the MHD equation of state was 
relevant only to the pre-main-sequence evolution. The MHD equation of state tables computed for this work 
considered H and He, plus a 13-element subset of the Grevesse & Noels (1993) heavy element composition 
(i.e., C, N, O, Ne, Na, Mg, Al, Si, S, Ar, K, Ca, and Fe); besides the effects due to neutrals and ions, effects 
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due to H _ , H2, and H^ were accurately taken into account, and approximate effects of the molecules C2, N2, 
O2, CH, CN, CO, NH, NO, OH, CO2, and H2O (although these molecules have no effect for solar models: 
they were actually added in anticipation of use in asymptotic giant branch models) . The OPAL equation of 
state considers only H, He, C, N, O, and Ne (plus hydrogen-molecule effects); however, the equation of state 
is quite insensitive to the precise makeup of the metallicity, and the OPAL equation of state may well be 
the more accurate one in the region where both are valid (Rogers et al. 1996; Gong, Dappen, & Nayfonov 
2001a). 

It has been noted that patching together equations of state can introduce spurious effects (Dappen et 
al. 1993; Basu, Dappen, & Nayfonov 1999). Such spurious effects were minimized in the present work by 
(1) performing a gradual switchover over a finite density or temperature region, to avoid any discontinuities, 
and (2) choosing switchover regions where the two equations of state were quite similar. The three switchover 
regions comprised —1.5 > log p > — 2 (corresponding to 5.8 > log T > 5.5 in the present Sun), 4.0 > log T > 
3.9, or 3.75 > log T > 3.7; in the first and last of these regions, the MHD and OPAL equations of state under 
solar conditions differ typically by several parts in 10 3 (with differences in the pressure P being several times 
smaller than this), i.e., only a few times worse than the expected precision of the equations of state. Near 
log T = 4, the differences were nearly an order of magnitude smaller (and also remained small over a very 
wide range of densities). Any spurious effects from the switchover should be smaller than spurious effects 
arising from the fact that the OPAL equation of state is not self-consistent, with the size of the inconsistency 
(in the outer Solar envelope) being comparable to the size of the differences between the OPAL and MHD 
equations of state (as discussed in § 3.1). 

Opacities: For interior temperatures (log T > 4), our reference standard solar model used the 1995 OPAL 
opacities (Iglesias & Rogers 1996) — these opacities use the Grevesse & Noels (1993) "GN93" solar com- 
position, and we refer to these opacities as "kopal:GN93-" The online opacity computation feature of the 
OPAL web page 3 allowed computation of OPAL opacities appropriate to the more recent Grevesse & Sauval 
(1998) "GS98" mixture ("kopal:GS98"); these were tested in variant models, as were opacities appropriate 
to the older Grevesse (1984) "Gr84" mixture ("«oPAL:Gr84")- The even older and much less accurate Los 
Alamos (LAOL) opacities "klaolss" (Keady 1985, private communication) were also tested in one case. At 
cool envelope temperatures (log T < 4), our reference standard solar model used the Alexander & Ferguson 
(1994) opacities "^Alexander" (which include molecular opacities). We also tested the effect of using the Sharp 
(1992) molecular opacities "fesharp" at cool envelope temperatures instead. 

In our reference standard solar model, we did our best to account for temporal and spatial variations in 
the opacity due to composition changes from diffusion and nuclear burning. The OPAL opacity tables allow 
interpolation of the opacity as a function of the hydrogen abundance and the metallicity Z; the abundances of 
the metals comprising Z are always proportional to Z in the OPAL tables, i.e., a "scaled solar" distribution. 
In addition, the OPAL opacity tables contain mixtures with excess carbon and oxygen (beyond that contained 
in the scaled solar metallicity), allowing interpolation in carbon and oxygen abundances. 

To take into account the variations in abundance due to diffusion, the metallicity value Z K that we 
used for metallicity interpolation in the OPAL opacity tables was proportional to the abundances of the 
elements heavier than oxygen: Z K = Z h , where Z h = Z a [J2 heavy X heavy i x i)o], where Z a and (X;) are 
the protosolar metallicity and composition, respectively, and ^2 heavy refers to a sum over elements heavier 
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than oxygen. In other words, we scaled the initial solar metallicity by the shift in heavy element abundances 
resulting from diffusion (note that our diffusion routines assumed that all metals diffused alike). However, for 
the CNO elements, there are additional changes due to nuclear burning, so that the CNO abundance profiles 
are not proportional to the heavy element abundance profiles; Turcotte et al. (1998) find that conversion of 
C and O into N results in an opacity change of ~ 1% that cannot be modelled by a change in Z alone. These 
variations in the CNO elements relative to Zh were accounted for in an approximate manner by an additional 
two-dimensional interpolation in nominal "excess carbon and oxygen" abundances C ex and O ex , where these 
excess abundances account for any variation in the CNO abundances relative to the scaled solar metallicity Z K 
of the OPAL opacity tables. Since there were no explicit opacity tables for changes in nitrogen, the best 
one could do was to distribute excess nitrogen equally between C ex and O ex (i.e., to assume that nitrogen 
opacities were midway between those of carbon and oxygen): C ex = C — C Z K /Z + 0.5(N — N Z K /Z ) and 
O ex = O - O ZJZ + 0.5(N - N Z K /Z ), such that CO e , = C ex + O ex = Z- Z K . 

Note that in general, as first carbon and then oxygen is burned to nitrogen, either C ex or O ex will 
be negative, implying extrapolation of the OPAL tables in the direction of zero C or O by a non-negligible 
fraction of Z. One might consider the uncertainties in such an extrapolation to be worse than the error 
inherent in treating all CNO opacities alike; in this case, if one of C ex or O ex was negative, one would set it 
to zero and subtract an equivalent amount from the other (so that the sum C ex + O ex remains unchanged) 
— this was the method followed in our reference models. However, we also tested the case where negative 
values were allowed, extrapolating as required. These turned out to yield essentially identical results to the 
reference model: the rms difference was only 0.00006 in the sound speed profile Ac/c and 0.0003 in Ap/p, 
comparable to the estimated numerical accuracy of the models. An alternative case where this latter CNO 
variation was not approximated by excess CO at all, i.e., having Z K = Zh but GQ ex = 0.0, also yielded 
almost identical results (rms difference of 0.00005 in Ac/c and 0.0004 in Ap/p). Similarly negligible rms 
differences were found for several cases testing different prescriptions for defining Zh, C ex , and O ex (differing 
in whether mass fractions or number densities were used, in how the excess nitrogen was divided between 
C ex and O ex , and in whether negative values were allowed for C ex and O ex ). 

Finally, one run was computed where CNO-interpolation in the opacities was actually performed. In 
addition to the standard OPAL opacity table, the OPAL web page 4 was used to compute an opacity table with 
a composition where C had been converted to N, and one where both C and O had been converted to N; the 
model CNO abundances were used to interpolate in opacity among these tables. This "CNO-interpolation" 
model also yielded results essentially identical to the reference model (an rms difference of 0.00006 in Ac/c 
and 0.0005 in Ap/p). 

An estimate for the upper limit of the effects of opacity uncertainties was made by simply setting 
Z K = f Z for a constant factor / = 0.9, 0.95, 1.0, 1.05, or 1.10; note that the case / = 1.0 is only slightly 
different from the reference standard solar model. 

Nuclear Reaction Rates: The reference standard solar model used the NACRE nuclear reaction rates 
(Angulo et al. 1999), supplemented by the 7 Be electron capture rates of Gruzinov & Bahcall (1997) for 
log T > 6 and of Bahcall & Moeller (1969) for log T < 6 (note that this latter low-temperature region is 
irrelevant for solar models). Variant models tested cases with nuclear rates changed according to the upper 
or lower limits quoted in the NACRE compilation; one case tested the use of the older Caughlan & Fowler 
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(1988) nuclear rates. The program uses the minimum of the weak (Salpeter 1955), intermediate (Graboske 
ct al. 1973), or strong (Itoh et al. 1979; Ichimaru & Utsumi 1983) screening factors (note that both weak and 
intermediate screening formulae include a term to take partial electron degeneracy into account) — for solar 
conditions, this choice leads to the use of weak screening, which is a very good approximation to the exact 
quantum mechanical solution there (see, e.g., Bahcall, Chen, & Kamionkowski 1998b; Gruzinov & Bahcall 
1998). Deuterium was not considered separately (it was assumed to have been entirely burned to 3 He on the 
early pre-main-sequence) , but all the other 15 stable isotopes up to and including 18 O were considered in 
detail (i.e., nuclear equilibrium was not assumed for any of them). The other stable isotopes up to 28 Si were 
included in the code (plus a few long-lived unstable isotopes, as well as Fe and a category for the sum of the 
other elements heavier than Si), but their nuclear reactions were not included since there are no significant 
effects under solar conditions (except for 19 F, which was assumed to be in CNO-cycle nuclear equilibrium for 
nuclear rate purposes). Neutrino capture cross sections were taken from Bahcall & Ulrich (1988), except for 
the 8 B-ncutrino cross section for capture on 37 C1, where the more recent value (5% higher) of Aufdcrheide 
et al. (1994) was used. 

Diffusion: A set of subroutines 5 were kindly provided to us (M. H. Pinsonneault 1999, private commu- 
nication) that take into account the diffusion (gravitational settling) of helium and heavy elements relative 
to hydrogen (see also Thoul, Bahcall, & Loeb 1994; Bahcall, Pinsonneault, & Wasserburg 1995). These 
subroutines assume that all heavy elements diffuse at the same rate as fully-ionized Fe; this yields surpris- 
ingly accurate results, as may be seen by comparing to the results of a more detailed treatment (Turcotte 
& Christenscn-Dalsgaard 1998; Turcotte et al. 1998). The upper limit of the effects of uncertainties in the 
diffusion constants was estimated by simply multiplying the diffusion constants, either for helium or for the 
heavy elements, by a constant factor. 

Convection: The Schwarzschild criterion was used to define convective boundaries; no core overshooting 
or envelope undershooting was allowed. Note that Morel, Provost, & Berthomieu (1997) found that including 
convective core overshooting had a negligible effect on the solar sound speed and density profiles, but on 
the other hand that including convective envelope undershooting by even a tenth of a pressure scale height 
moved the solar convective envelope boundary inwards by eight times the uncertainty in the observed value, 
yielding a sharp spike in the difference between observed and calculated sound speed profiles. Rotation- 
induced mixing was not considered in our models; the effect that it would have is discussed in § 3.1. 

Composition: The reference standard solar model used a value of Z/X = 0.0245 for the the present solar 
surface mctals-to-hydrogen ratio (by mass fraction), as given by Grevesse & Noels (1993). Variant models 
tested the ~ 13% higher older value of Z/X = 0.0277 (Grevesse 1984; Anders & Grevesse 1989), as well as 
the 6% lower value of Z/X = 0.023 recommended as a "preliminary" value by the recent work of Grevesse 
& Sauval (1998). Values of Z/X = 0.0203 and 0.0257 were also used, for cases where the C, N, O, and Ne 
abundances of the Grevesse & Sauval (1998) mixture were all either decreased or increased, respectively, by 
their quoted uncertainties. 

For a given solar model, the presolar abundances of the heavy elements were always taken from the 
OPAL opacity tables that we used for that solar model. For the reference standard solar model and most 
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variants, the "kopal:GN93" opacities were used, which have the composition mix of Grevesse & Noels (1993). 
A variant model with a low value Z/X = 0.023 used instead the composition mix of Grevesse & Sauval (1998), 
and the corresponding "kopal:GS9s" opacities. Similarly, a variant model with a high value of Z/X = 0.0277 
used of the abundance mix of Grevesse (1984), and the corresponding "KoPAL:Gr84" opacities. Tests were 
also made with C, N, O, and Ne abundances increased or decreased by their uncertainties of 15% relative 
to Fe, with the OPAL opacities appropriate to these revised mixes: we refer to these as "kopal:GN93TC-Nc" 
and "kopal:GN931C-Nc," respectively, when variations are relative to the Grevesse & Noels (1993) mix, and 
as "kopal:GS98TC-Nc" and "kopal:GS981C-Nc" when variations are relative to the Grevesse & Sauval (1998) 
mix. 

Solar Mass: A present solar mass of M© — 1.9891 x 10 33 g (Cohen & Taylor 1986) was used in all cases. 
Their quoted uncertainty of 0.02% is too small to have any significant effect, and is in fact smaller than the 
amount of mass lost by the Sun in the form of radiation (AM = AE/c 2 yields a mass loss of 0.03%, where 
AE is the total energy radiated away via photons and neutrinos since the Sun formed). Mass loss from the 
present solar wind is also small, and even if the average solar wind over the past 4.6 Gyr had been an order 
of magnitude higher than its present value, as is suggested by measurements of noble gas isotopes in lunar 
rocks (Geiss 1973; Geiss & Bochsler 1991; Kerridge et al. 1991), the total amount of mass lost would be 
about 0.2% (i.e., a total mass loss of < 0.002 M© during the Sun's lifetime up to the present). Mass loss of 
this amount would yield negligible changes in the solar sound speed profile (about a part in 10 4 ), as shown 
in our companion paper "Our Sun V" (Sackmann & Boothroyd 2002). Solar mass loss was therefore ignored 
for all cases considered in this paper. 

Note that, based on a correlation of X-ray flux with their measured mass loss rates in nine GK dwarfs, 
Wood et al. (2002) have recently proposed a mass loss time dependence M oc t~ 2 00±0 - 52 in such stars (i.e., 
in main sequence stars with masses not too far from that of the Sun), with a maximum mass loss rate 
of ~ 10 3 times that of the present solar wind. This is discussed in our companion paper "Our Sun V" 
(Sackmann & Boothroyd 2002), and will be investigated in more detail in a future work (A. I. Boothroyd & 
I. -J. Sackmann, in preparation). In summary, total solar mass loss from the formula of Wood et al. (2002) 
could be of order 0.01 M©, but most of this would take place very early on the main sequence due to the 
fr 2 time dependence, so the effect should be relatively small, at most a few parts in 10 4 in the sound speed 
profile. 

Solar Luminosity: For the reference standard solar model and most variant cases, a present solar lumi- 
nosity of L© = 3.8 54 x 10 33 erg s _1 was used, as discussed in Sackmann et al. (1993). A value 0.3% lower 
(3.842 x 10 33 erg s _1 ) with an estimated 1-a uncertainty of 0.4% was recently obtained by Bahcall et al. 
(2001), based on the observations of Frohlich & Lean (1998) and Crommelynck et al. (1996). Variant models 
considered the effect of using this more recent solar luminosity value, and high and low values 2- a (0.8%) 
above and below it. Note that most of the uncertainty in Lq comes not from uncertainties in the present 
solar irradiance, but rather from uncertainties in the slight long-term variability of the solar luminosity. 
For example, Lean (2000) estimates a difference of 0.2% between the present value and that during the 
seventeenth century Maunder Minimum. 

Solar Radius: The reference standard solar model used a solar radius at the photosphere (r = 2/3) of 
Rq = 695.98 Mm (Ulrich & Rhodes 1983; Guenther et al. 1992). Variant models considered the effect of 
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using the value of 695.78 Mm suggested by the helioseismic /-mode study of Antia (1998), or the value of 
695.508 Mm suggested by the solar-meridian transit study of Brown & Christensen-Dalsgaard (1998). 

Solar Age: In his Appendix to Bahcall et al. (1995), G. J. Wasserburg provides a systematic analysis of 
the upper and lower bounds on the age of the Sun, as obtained from isotopic ratios measured in meteorites. 
We briefly paraphrase his discussion in this paragraph: The protosolar nebula (out of which the Sun and the 
meteorites formed) contained not only the stable isotope 27 A1 but also the unstable isotope 26 Al, which decays 
into 26 Mg with a half-life of only 0.7 Myr. This 26 Al must have been injected into the protosolar nebula 
from the stellar source where it was created. Isotopic measurements of mctcoritic crystallized refractory 
condensates show that they had a ratio 26 A1/ 27 A1 = 5 x 10~ 5 at the time they formed. Even if the stellar 
source of the 26 Al had a very high ratio 26 A1/ 27 A1 ~ 1, the time interval between the formation of the 26 Al 
and the formation of the meteorite cannot have exceeded ^11 Myr (due to the short decay timescale of 26 Al). 
The Sun cannot have formed later than these meteorites. However, the Sun/meteorite system must have 
formed after the injection of 26 Al into the protosolar nebula, so the Sun cannot have formed earlier than 
^11 Myr before the formation of these meteorites. The age of the meteorites has been accurately measured 
using 207 Pb/ 206 Pb ratios, to be 4.565 ± 0.005 Gyr. It follows that the Sun cannot have formed earlier than 
4.59 Gyr ago, nor later than 4.55 Gyr ago (i.e., 4.565 + 3 x 0.005 + 0.011 Gyr, or 4.565 - 3 x 0.005 Gyr, at 
the 3-cr level). This total solar age estimate of, in effect, t Q — 4.57 ±0.01 Gyr is in agreement with the total 
solar (and solar system) age t ss = 4.53 ± 0.03 Gyr inferred previously by Guenther (1989) from a formation 
age of 4.53 ± 0.02 Gyr (Wasson 1985) for the meteorites and planets. 

Our solar models were started relatively high on the pre-main-sequence Hayashi track, with central 
temperatures below 10 6 K; note that the ages of all our models are quoted relative to this pre-main- 
sequence starting point, which is within a few Myr of the solar formation age constrained by the meteoritic 
ages discussed above. A relatively high value of i Q = 4.6 Gyr was used for the reference standard solar 
model (cf. the observationally inferred value t Q — 4.57 ± 0.01 Gyr from the previous paragraph). To get a 
reliable estimate of the sensitivity of the models to the solar age uncertainty, variant models were computed 
with ages differing by very large amounts, namely, i = 4.5 and 4.7 Gyr. 

From our Hayashi track starting point, it takes only ~ 3 Myr for the luminosity on the pre-main-sequence 
to drop below 1 L & (i.e., below the present solar luminosity), but much longer, namely an additional 40 Myr, 
to reach the zero-age main sequence (ZAMS) — we have defined the ZAMS as the stage where the pre-main- 
sequence contraction terminates and the Sun begins to expand slowly, as nuclear burning in the core (rather 
than gravitational contraction) supplies essentially all of the solar luminosity. For the next ~ 50 Myr on 
the early main sequence, evolution is fairly fast, as p + C reactions burn the initial carbon to nitrogen near 
the Sun's center, resulting in a short-lived convective core. Subsequently, the central carbon and nitrogen 
abundances approach their CN-cycle equilibrium values, the convective core disappears, and the Sun settles 
down to burn hydrogen mainly via the pp-chain reactions. Note that the pre-main-sequence timescale 
implies that the total solar age t Q used in this paper can be converted into a main sequence solar lifetime 
by subtracting about 0.04 Gyr — this was also pointed out by Guenther (1989). 

Helioseismology: We compared our solar models to profiles of the solar sound speed cq, density p Q , and 
adiabatic index (Fi)q obtained from the helioseismic reference model of Basu et al. (2000) 6 , which they 



6 From the denser-grid machine-readable form of their Table 2, at http://www.sns.ias.edu/~jnb/ 
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obtained by inversion from the helioseismic frequency observations. In the inversion process, a standard 
solar model is required, but Basu et al. (2000) demonstrated that the resulting c© and p© profiles of the 
helioseismic reference model are relatively insensitive to uncertainties in the standard solar model used for 
this purpose (except for uncertainties in Rq, as discussed in § 3.1). They estimated a net uncertainty of few 
parts in 10 4 for the sound speed c© and adiabatic index (Fl)©, and a few parts in 10 3 for the density p©. 
However, in the Sun's core (r < 0.1 Rq), systematic uncertainties in the helioseismic sound profile are 
increased by a factor of ~ 5; this was demonstrated by Bahcall ct al. (2001), who compared helioseismic 
inversions of different helioseismic data sets. We used their comparison to estimate the r-dependence of 
the systematic error in c© in the core and in the convective envelope (namely, a fractional systematic error 
decreasing linearly from 0.0013 at r = 0.05 Rq to 0.0003 at r = 0.2 Rq, constant from there to r = 0.72; Rq, 
then increasing linearly to 0.00052 at r = 0.94; Rq). For c©, this systematic error can be significantly larger 
than the statistical errors quoted in the Table 2 of Basu et al. (2000), and we combined the two in quadrature 
to get the fractional error (<r c /c) for the purpose of calculating weighted rms differences — the rms fractional 

difference in c is given by [(^ c / c )/( fT c/c)] 2 | / [l/(cr c /c)] 2 |^ . For (Ii)© and p©, the systematic 

errors are comparable to or smaller than the statistical ones, and the statistical errors sufficed for calculating 
weighted rms differences. 



3. Results and Discussion 

3.1. Sound Speed and Density Profiles 

We present in Figures 1 through 7 profiles of the adiabatic sound speed differences 5c/ c = (c — 
c modei) I ' c ©; profiles of the density differences Sp/p = (p© — p mo dei) I Pq are available online 7 . For our equation 
of state comparison, we also considered the equivalent fractional difference in the adiabatic index Ti . Note 
that we use "<5" to denote differences between the helioseismic profile and one of our models, and "A" to 
denote differences between two of our models with different input parameters — the "<5" values are the 
profiles plotted in our figures, while the "A" values refer to the difference between one plotted curve and 
another. 

The theoretical sound speeds c mo dei and densities pmodei are from our computed reference standard 
solar model and from our variant standard solar models. Our reference standard solar model used current 
input parameters, as discussed in § 2; our variant standard solar models comprised standard solar models 
with one or more input parameters varied within the permitted range. We present all our sound speed and 
density profiles in terms of differences relative to the observed helioseismic reference profiles of Basu et al. 
(2000). This choice of presentation not only allows one to see the effects of the uncertainties in the input 
parameters, but also shows which choice of input parameters agrees best with the helioseismic observations. 
Our "OPALeos-lowT" model (discussed in more detail below) is the one most nearly comparable to the 
"STD" model of Basu et al. (2000), and yields similar 5c/ c and <5ri/ri curves, as may be seem by comparing 
their Figures 2 and 3 to our Figure 1. Their models fit the solar sound speed profile somewhat better than 
our models do, but not significantly so, considering the size of the effects (discussed in detail further below) 
that result from reasonable variations in the input parameters of the solar model. We made no attempt to 
compare our solar models in the region outside r = 0.943 Rq, the last point on the helioseismic profiles of 
Basu et al. (2000) — the reason that this is their outermost point is that significant systematic uncertainties 



7 http://www. krl.caltech.edu/~aib/papdat. html 
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arise in inversions near the solar surface (see, e.g., Di Mauro et al. 2002). However, we have compared the 
OPAL and MHD equation of state in this region, as have some other investigators (see, e.g., Guzik & Swenson 
1997; Richard et al. 1998; Basu et al. 1999; Gong et al. 2001a; Di Mauro et al. 2002); this is discussed in 
detail further below. 

Rotation effects: The prominent peak in Sc/c at r ~ 0.7 Rq visible in Figure la is due to the neglect 
of rotation-induced mixing just below the base of the solar convective envelope, as has been shown by a 
number of investigators who have included parameterized rotational mixing (see, e.g., Richard et al. 1996; 
Brun, Turck-Chieze, & Zahn 1999; Basu et al. 2000; Bahcall et al. 2001; Turck-Chieze et al. 2001b). There 
are significant uncertainties in the physical processes that lead to rotation-induced mixing. However, all 
of these investigators agree that rotational mixing is capable of smoothing out the peak at r ~ 0.7 Rq, 
and that this has a relatively small effect on the sound speed elsewhere in the Sun (a fractional change 
of less than 0.001). For example, Bahcall et al. (2001) found that including "maximal" rotational mixing 
spread out this peak over the region 0.3 Rq < r < 0.7 Rq, eliminating the prominent peak but worsening 
the agreement with the helioseismological sound speed profile by about 0.001 in much of the solar interior 
(0.3 Rq < r < 0.6 Rq). The "minimal" mixing model of Richard et al. (1996) yielded much the same result, 
as did similar models of other authors (Gabriel 1997; Brun et al. 1999; Turck-Chieze et al. 2001b). As far 
as the sound speed and density profiles in the core and the convective envelope are concerned, rotational 
mixing below the base of the convective envelope should have no significantly effect, as shown by the above 
authors. 

Since the prominent peak at r ~ 0.7 Rq results from the neglect of rotational mixing, we did not 
require agreement in this region between profiles from our theoretical models and profiles inferred from 
the helioseismic observations. Nor did we require agreement in the inner core region, since the present 
helioseismic observations still result in large uncertainties in the inferred profiles there; for example, as 
shown by Bahcall et al. (2001), the use of a different helioseismic dataset could remove the the sharp upturn 
in Sc/c for r < 0.1 Rq in Figure la (and even convert it into a downward trend). On the other hand, we 
aimed for agreement in the regions 0.1 Rq < r < 0.6 Rq and 0.72 Rq < r < 0.94 Rq, where disagreements 
are due to imperfections in the input physics or uncertainties in the observed solar parameters. This is 
demonstrated by our variant models, and by the variant models of other investigators (see, e.g., Morel et al. 
1997; Guzik & Swenson 1997; Richard et al. 1998; Basu et al. 2000; Gong et al. 2001a; Neuforge-Verheecke 
et al. 2001a,b). 

Convergence accuracy effects: The accuracy with which the model is converged to the solar radius, 
luminosity, and Z/X values can affect the sound speed profiles. The extent to which this occurs depends 
on the accuracy of the convergence, and on the sensitivity of the profiles to variations in Rq, Lq, and Z/X 
(these are discussed in detail below) . Our convergence accuracy resulted in effects no larger than a few parts 
in 10 5 on the sound speed in the solar interior (r < 0.6 Rq). In the convective envelope, where the sound 
speed is quite sensitive to Rq, the effect was typically less than a part in 10 4 , but could be as high as a few 
parts in 10 4 for the few models with the worst convergence in Rq. 

Note that the sound speed profile in most of the solar convective envelope (0.72 Rq < r < 0.94 Rq) 
is sensitive mainly to the equation of state and to the solar radius, with other uncertainties having only a 
minor effect there, as can be seen by considering this subregion in Figures 1 through 7. 



- 11 - 




-0.003 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

radius r (R ) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
radius r (R@) 



Fig. 1. — Changing the zoning, the outer-envelope equation of state, or the solar luminosity: the effects 
on (a) the adiabatic sound speed c, and (b) the adiabatic index I\. The reference standard solar model 
(thick solid line: errorbars give statistical error only, for the inferred helioseismic profile) switches from the 
OPAL to the MHD equation of state for log p < -2 (i.e., r > 0.94 R Q or log T < 5.5). The "OPALcos- 
midT" (short-dashed line) and "OPALeos-lowT" (dashed line) models have the equation-of-state switchover 
at log T « 4.0 and log T w 3.75, respectively. The thin solid lines show the effects of using the maximum 
("Lhigh") and minimum ( u Li ow ") values of L©; for clarity, these are shown only for the region r < 0.3 Rq 
where there is a visible effect. 



- 12 - 



Zoning effects: We investigated the effects of using two different zonings. Our coarse-zoned models had 
about 2000 spatial zones in the model, and about 200 time steps in the evolution from the zero-age main 
sequence to the present solar age (plus about 800 time steps on the pre-main-sequence) . These values are 
comparable to those used by most authors for solar models, although many authors use more complex 
algorithms than we did to compute changes between one timestep and the next (allowing the use of fewer 
main-sequence timesteps at the cost of more CPU-time per timestep), and some authors ignore the pre- 
main-sequence evolution entirely (or use less stringent accuracy conditions there). Typically, we ran solar 
evolutionary sequences (iteratively improving the input parameters Z , Y , and a) until our coarse-zoned 
models were converged to match the solar luminosity and radius to about a part in 10 5 , and the solar 
surface Z/X to a part in 10 ; a few cases where convergence was slow were nearly 10 times worse. Our fine- 
zoned models had 10 000 spatial zones and took 1500 main-sequence time steps (plus 6000 pre-main-sequence 
time steps) — a factor of 5 increase in both spatial and temporal precision — and were typically converged 
to better than a part in 10 5 for Rq and L Q , and a few parts in 10 5 for Z/X. We also tested some very 
coarse-zoned models, with 1000 spatial zones, 100 main-sequence time steps (plus 600 pre-main-sequence 
time steps), and convergence to the solar parameters of a few parts in 10 4 . (A coarse-zoned converged solar 
model took a few hours of CPU-time on a fairly high-performance ES40 computer, as compared to a few days 
of CPU-time for a fine-zoned converged model; these times were roughly tripled on a 450 Mhz Pentium III 
PC.) 

Figure la shows that the fine zoning made only a very modest improvement relative to the coarse-zoned 
case, less even than the statistical errors in the sound speed and density profiles obtained from helioseismic 
inversions. Even the very coarse-zoned test case did not do too badly: in the solar interior, it differs from the 
fine-zoned case by no more than 0.0004 in the sound speed profile and 0.003 in the density profile (amounts 
comparable to the systematic uncertainties in the helioseismic inversion) — the rms differences (over the 
entire Sun) are even smaller, namely rms{Ac/c} w 0.0003 and rms{Ap/p} « 0.002. The coarse-zoned model 
did about 3 times better still, with rms differences relative to the fine-zoned case of rms{Ac/c} « 0.0001 
and rms{Ap/p} « 0.0008. Zoning changes had no effect on the adiabatic index Pi (the coarse-zoned Ti 
curves were not plotted in Fig. lb, since they would be precisely superimposed on the fine-zoned curves). 
Additional tests demonstrated that changes in the coarseness of zoning always led to the same negligibly 
small systematic shift in solar interior sound speed and density values (although inaccuracies in matching the 
observed solar surface parameters could lead to slightly larger random variations in the convective envelope 
region r > 0.7 Rq). We therefore felt justified in running most of the models with our coarse zoning. Note 
that Morel et al. (1997), with about 1000 spatial zones, 60 main-sequence time steps, and convergence to 
present solar surface parameters of a part in 10 4 (similar to our very-coarse-zoned case), claimed a numerical 
internal accuracy of 0.0005 in the sound speed, similar to what we found for our very-coarse-zoned case. 

Equation-of-state effects: Gong et al. (2001a) compared four current equations of state: their own 
MHD equation of state (Dappen et al. 1988), the OPAL equation of state (Rogers et al. 1996), the CEFF 
equation of state (Christensen-Dalsgaard & Dappen 1992), and the SIREFF equation of state Guzik & 
Swenson (1997). In each case, they used the approximation of a 6-element composition mixture (H, He, 
C, N, O, Ne) with the same composition as in the OPAL equation of state of Rogers et al. (1996). For 
3.7 < log T < 6, at corresponding solar densities but without the usual "r-correction" that eliminates the 
short-range divergence in the Debye-Huckel potential, Gong et al. (2001a) find OPAL — MHD differences in 
the equation of state of AP/P < 0.0005, A Xp < 0.005, A X t < 0.007, and APi < 0.004, with differences 
several times smaller at 6 < log T < 7, and comparable differences between other pairs of equations of state 
- recall that Xp = {d\n P/d\n p)t, Xt = (91nP/<91nT) p , and the adiabatic index is Ti = (d\nP/dh\p) s . 
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With the usual "T-correction," but for a pure hydrogen- helium mixture (Z = 0), they find much larger 
differences: AP/P < 0.002, A Xp < 0.02, A X t < 0.03, and AI\ < 0.008. 

Guzik & Swcnson (1997) compared the effect in different solar models, finding roughly the same OPAL — 
MHD difference as the last of the above comparisons (namely, the case with the "r-correction" ) , except that 
they find a much larger difference in the pressure; in the region r > 0.95 Rq (i.e., log T < 5.5), they find 
AP/P < 0.014, AC p /C p < 0.036, and AI\ < 0.006. Note that these differences include the effects of slightly 
different temperature, density, and composition profiles in the different solar models (which may cither 
increase or decrease the differences in the thermodynamic quantities, since the solar models are re-adjusted 
to reproduce the present solar luminosity, radius, and surface composition). 

We performed our own OPAL — MHD comparison, at fixed temperature, density, and composition grid- 
points in the equation of state tables — the same comparison as that performed by Gong et al. (2001a). Our 
results are consistent with those of Guzik & Swenson (1997); for a typical solar {T, p} profile in the outer 
envelope region, we found AP/P < 0.017, AC V /C V < 0.03, A Xp < 0.015, A X t < 0.032, and Ari/Ti < 0.006. 

The MHD equation of state is obtained in a fully self-consistent manner from the free energy; inaccuracies 
can arise only from deficiencies in the formulas used to obtain the free energy (Dappen et al. 1988; Gong, 
Dappcn, & Zedja 2001b). However, for the OPAL equation of state (Rogers et al. 1996), we found that 
there were significant inconsistencies when we compared their tabulated values of Ti, T2HT2 — 1), and 
(r 3 — 1) to values calculated from their tabulated values of P, C v , Xp , and \t- In the solar core, these 
inconsistencies are very small (a few parts in 10 4 ), but in the outer envelope (r > 0.94 Rq, or log T < 5.5) 
we found inconsistencies as large as 3% at grid-points that would be used in the OPAL interpolation formulae 
when computing thermodynamic quantities (although the grid-points nearest to the solar {T, p} locus have 
inconsistencies of less than 1%). Generally, the size of these inconsistencies varied smoothly in the OPAL 
grid, but in at least a few positions a few grid spacings away from the solar {T, p} locus, there were sharp 
"spikes" where one of the thermodynamic quantities had a "glitch" (an error of several percent) at a just a 
couple of adjacent density and/or temperature points. In addition, for 4 of the 8 lowest OPAL T-grid points 
(at 3.71 < log T < 3.77), there is a "sawtooth" error: at every second density value, the tabulated quantities 
are shifted systematically relative to the values at neighboring densities and temperatures. These shifts can 
be as large as 1% at low densities for C v and X t, and are always of order 0.1% for P. 

A new OPAL 2001 equation of state has recently become available 8 (see also Rogers 2000, 2001), which 
includes relativistic electron effects and extends to both lower temperatures and higher densities then the 
original OPAL equation of state. Preliminary tests indicate that this OPAL 2001 equation of state has larger 
but smoother inconsistencies in its tabulated thermodynamic quantities, with few or no glitches, except in 
the extended high-density region (where there are some very large ones) and at the 2 lowest T-grid points 
(log T < 3.35). 

Gong et al. (2001a) found that the changing from a 6-element composition mixture to a 15-element 
mixture has an effect an order of magnitude smaller than the OPAL — MHD differences (we also performed 
such tests, and came to the same conclusion). They nonetheless recommended the use of at least a 10-element 
mixture for the best accuracy in an equation of state. They also noted a couple of minor deficiencies in the 
MHD equation of state that "moves it away from both hclioseismically determined values and OPAL" (Gong 
et al. 2001a). 



ftp://www-phys.llnl.gov/pub/opal/cos2001/ 



-14- 



Morel et al. (1997) compared solar models with the OPAL equation of state to ones with the CEFF 
equation of state, and found a small but non-negligible effect: sound speed differences Ac/c of slightly over 
0.001 and density differences Ap/p up to 0.01. Note that they had set the value of Z eos used in their equation 
of state to a fixed value of 0.019; however, since the equation of state is only weakly sensitive to Z, this should 
have only a minor effect on their models, and should not affect their comparison of the two equations of state. 
In general, the temperature, density, and composition profiles would all be slightly different between solar 
models with different equations of state, since the input parameters Zq, Yq, and a are adjusted individually 
for each solar model to obtain the best fit to the present solar luminosity, radius, and surface composition. 

Guzik & Swenson (1997) presented more extensive solar model comparisons, comparing both their 
own SIREFF equation of state and the MHD equation of state to the OPAL equation of state (they too 
used a fixed Z eos , of 0.02, but again this should not affect the comparison). They likewise found an effect 
Ac/c < 0.001 at r < 0.95 Rq (with differences up to 0.004 near the surface). They also presented differences 
between the values of the pressure P, specific heat at constant pressure C p , internal energy U, and adiabatic 
index L\. As mentioned above, near the solar surface (r > 0.95 Rq) differences between models with 
different equations of state were relatively large, of order 1%. However, for r < 0.9 Rq Guzik & Swenson 
(1997) reported OPAL - MHD differences between their solar models of AP/P < 0.0015, AU/U < 0.002, 
Ali/ri < 0.0007, and AC P /C P < 0.006; the OPAL- SIREFF differences were slightly smaller for C p , slightly 
larger for P and U, and much larger (a factor of ~ 3) for L\. For r > 0.3 Rq, the SIREFF value of Fi has 
several relatively large "wiggles" (fractional variations ~ ±0.002) relative to either OPAL or MHD. On the 
other hand, most of the difference in the core, and perhaps some difference in the average trend further out, 
may be due to the fact that SIREFF includes relativistic electron effects, while MHD and OPAL do not — 
although they are included in the new OPAL 2001 equation of state (Rogers 2000, 2001), and have recently 
been added to the MHD equation of state by Gong et al. (2001b). 

Elliot & Kosovichev (1998) estimated that inclusion of relativistic effects would reduce the MHD or 
OPAL value of Ti by a fraction 0.002 at r w 0.1 Rq, this correction growing smaller with increasing r, to 
reach 0.001 at r « 0.3 Rq and zero near the solar surface. They pointed out that such a shift in L\ for 
models using the OPAL or MHD equation of state would significantly improve the agreement in the solar 
interior with the inferred helioscismic L\ profile. Certainly, if such a correction were applied to our L\ curve 
in Figure lb, the model profile would agree with the hclioseismic profile within the statistical errors for all 
r < 0.6 Rq (recall that decreasing a model quantity shifts the curve upwards in the figures). Recently, Gong 
et al. (2001b) confirmed that adding relativistic electron effects to the MHD equation of state yields a change 
in L\ very close to that estimated by Elliot & Kosovichev (1998). 

The adiabatic sound speed is defined as c = (FiP/p) 1 / 2 ; changes in the solar ratio of P/p would result 
not only from changes in the equation of state but also from readjustments of the solar structure in response 
to these changes, so it is not obvious a priori what effect relativistic corrections would have on the sound 
speed. Consideration of the effect from T\ alone suggest that relativistic corrections might reduce the slope 
at r < 0.5 Rq in the 5c/ c curve of Figure la. The sound speed differences presented by Guzik & Swenson 
(1997) for their OPAL — SIREFF comparison suggest that this would in fact be the case, and that a fractional 
decrease of order 0.001 in the sound speed c mo dei near the Sun's center (i.e., an increase of 0.001 in 5c/ c 
there) would result from relativistic corrections to the OPAL equation of state. 

Richard et al. (1998) compared solar models with the OPAL and MHD equations of state, looking at 
the value of Y\ in the convective envelope (0.72 Rq < r < 0.98 Rq); they found that the OPAL equation of 
state appeared to perform slightly better there. This is not very surprising; what is perhaps more surprising 
is how well the MHD equation of state does in the solar interior, since it was only originally designed to be 
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accurate for p < 10~ 2 g cm~ 3 (D. Mihalas 1999, private communication; see also Dappcn ct al. 1988; Gong 
et al. 2001a) — note that this p condition corresponds to r > 0.94 Rq and log T < 5.5 in the Sun. We 
investigated the effect of changing the equation of state only in this outer region where both are expected 
to be valid. Note that, while the OPAL opacity tables are unreliable for log T < 4 due to their neglect of 
molecular opacities, the OPAL equation of state tables include molecular hydrogen effects, and should be 
reasonably accurate down to their lower tabulation limit of log T — 3.699 (Rogers et al. 1996). 

Basu et al. (1999) considered the outer part of the solar convective envelope (0.8 Rq < r < 0.99 Rq); 
they used a hclioseismic inversion to study intrinsic Li differences relative to the solar values, where 
(STi/Ti)i nt is only that part of the difference that is ascribed to the equation of state in the inversion. 
They found that the OPAL L\ value was preferable to the MHD value for r < 0.97 Rq. Both equations of 
state did relatively well in the inner envelope, with (STi/Ti)i n t < 0.0005; but near the solar surface both 
yielded a relatively large (negative) peak with (dTi/Ti)i nt ~ 0.004. This peak was deeper in for the MHD 
equation of state (at r ~ 0.96 Rq, as compared to r ~ 0.975 Rq for OPAL), and also wider (reaching in 
to r ~ 0.9 Rq, as opposed to r ~ 0.95 Rq for OPAL). Di Mauro ct al. (2002) used recent observations 
of high-degree modes to perform a similar comparison, finding essentially the same result, though with a 
slightly higher peak difference (STi/Ti)i nt ~ 0.0055 near the surface. Our Figure lb plots the total L\ differ- 
ences rather than the intrinsic ones, but it is nonetheless likely that the sharp downturn in our MHD <5Fi/Fi 
profile (solid curve) at 0.9 Rq < r < 0.94 Rq corresponds to the inner edge of this "peak" in the MHD 
disagreement. 

Our reference standard solar model used the OPAL equation of state in the interior regions, switching 
over to the MHD equation of state in the outer envelope; this switchover was performed gradually over the 
region —1.5 > log p > — 2 (corresponding to 0.89 Rq < r < 0.94 Rq and 5.8 > log T > 5.5 in the present 
Sun) . We compared this reference standard model with two cases where the switchover occurred even further 
out in the envelope: a case "OPALeos-midT" where the switchover occurred for 4.0 > log T > 3.9, and a 
case "OPALeos-lowT" where the switchover occurred for 3.75 > log T > 3.7. In this latter "OPALeos-lowT" 
case, the MHD equation of state is used only outside the Sun's photosphere, and thus has negligible effect on 
the main sequence evolution. On the other hand, any artifacts induced by the switchover might be smaller in 
the former "OPALeos-midT" case, since differences between the two equations of state are significantly less 
near log T = 4 than near log T = 3.7 (although any such artifacts should be small in any case, comparable 
to effects of the inconsistencies in the OPAL equation of state, as discussed in § 2). 

Figure 1 shows that, as one would expect, changing the equation of state in the outer envelope alone 
has no effect on the interior, and in fact there is only a minor effect in that part of the convective envelope 
where the equation of state remains unchanged. (This also demonstrates that any artifacts from the equation 
of state switchover do not affect the solar interior, and also have at most a minor effect on the convective 
envelope.) Since the effects were so small, we computed fine-zoned cases for this equation-of-state test; these 
are the ones presented in Fig. 1. 

Note that L\ is affected significantly only by variations in the equation of state. 

Solar luminosity effects: Bahcall et al. (2001) tested the effect of 2-cr changes in the value of Lq (namely, 
±0.8%) on their solar models, finding only a minor effect on neutrino fluxes and negligible effects on the 
other quantities they considered. For completeness, we made the same test with our own models, confirming 
their results. Figure la shows the effects on the sound speed of a solar luminosity 0.8% lower ("Li ow ") and 
0.8% higher {"L high ") than the most recent value of 3.842 x 10 33 erg s' 1 (Bahcall et al. 2001; Frohlich & 
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Fig. 2. Effect on the sound speed of changing the solar radius or the solar age. The reference standard 
solar model (thick solid line) has the standard solar radius (Rq — 695.98 Mm: "R98") and a (high) solar 
age (i© = 4.6 Gyr). The dot-dashed line shows the effect of using a very low solar age t© = 4.5 Gyr (the 
observational value of t© = 4.57 ± 0.01 Gyr implies a total 3 — a age range only half as large as the difference 
between these cases). The thick long-dashed line shows the effect of using the low value of Rq = 695.78 Mm 
("R78"), while the thick short-dashed line shows the effect of using the very low value of Rq — 695.508 Mm 
( "R508" ) . The dotted lines show the same radius comparisons for the "OPALeos-midT" case (in which the 
MHD equation of state is used only for log T < 4.0, rather than log T < 5.5). Note that variations in the 
solar radius would also affect the inferred helioseismic sound speed profile; this is discussed in the text, but 
the effects are not included in this figure. 

Lean 1998; Crommelynck et al. 1996) — note that our reference standard solar model lies closer to Lhi g h 
than to Li ow , as it uses a slightly higher Lq value than the most recent estimate (see § 2). To avoid confusion 
with other curves, the "Lhigh" and "Lj ou) " curves are shown only in the region where they differ the most, 
namely, r < 0.3 Rq; even in this region, a shift of 0.8% in Lq produces a fractional change in the sound 
speed of less than 3 parts in 10 4 , dropping to 1 part in 10 4 for r > 0.3 Rq. 

Solar radius effects: Basu et al. (2000) demonstrated that using a solar radius different from the standard 
value of Rq — 695.98 Mm could have a small but not completely insignificant effect on both the sound speed 
profile inferred from helioseismic inversions and that computed in solar models. They found that using 
the 0.03% smaller solar radius value Rq = 695.78 Mm (case "R78") suggested by the /-mode study of 
Antia (1998) would reduce the inferred helioseismic sound speed profile throughout the Sun by about the 
same small fraction, namely 0.0003 (nearly independent of position in the Sun); similarly, using the 0.07% 
smaller value Rq = 695.508 Mm (case "R508") suggested by the solar-meridian transit study of Brown & 
Christensen-Dalsgaard (1998) would reduce the inferred helioseismic sound speed profile by 0.0007. On plots 
such as ours of fractional differences 5c/ c = (c© — c m odei) / cq 1 such a reduction in "c©" would shift all the 
curves downwards by the given amounts. This shift has not been performed in Figure 2 — we only display 



-17- 



our c mo dei values relative to the c Q values of Basu et al. (2000), so as to allow comparisons between different 
solar models — but this effect has been included in our quoted rms values relative to the Sun "rms{<5c/c}" 
(see next paragraph, and Table 1). Basu et al. (2000) also calculated that a change in the solar radius would 
result in a non- uniform shift in the inferred hclioseismic p© and (Ii)© profiles, by amounts comparable to 
the statistical errors in these quantities; although these shifts are barely significant statistically, in contrast 
to the shift in c©, for completeness their effects have been applied to our rms values calculated relative to 
the helioseismic profiles for the "R78" and "R508" cases. 

Figure 2 illustrates the effect on the sound speed of changing the value of Rq from the standard "R98" 
case to the smaller "R78" and "R508" cases. Only in the convective envelope is the sound speed significantly 
affected, with the largest effect being near the solar surface. In the "peak" region just below the convective 
envelope, these "R78" and "R508" curves differ from the reference standard solar model by no more than 
a few parts in 10 4 , and this difference drops to about a part in 10 5 for r < 0.6 Rq. In the convective 
envelope (r > 0.72 Rq), the "R78" case is an improvement on the reference model (reducing envelope-only 
rms env {Sc/c} from 0.0007 to 0.0004 when one includes the effect of the shift in the inferred helioseismic 
profiles), but the "R508" case is worse (rms em ,{<5c/c} = 0.0011). For the "OPALeos-midT" case, reducing 
the solar radius always worsens agreement in the convective envelope (rms en „ {Sc/c} of 0.0003 is increased to 
0.0009 or 0.0018 for "R78" or "R508," respectively). Figure 2 also shows that the effect on the sound speed 
profile of changes in the solar radius adds linearly to effects from changes in the envelope equation of state. 

The overall rms and the rms in the interior are not much affected by variations in the solar radius, as 
may be seen from Table 1. However, if relativistic corrections had been included in the equation of state, 
the Sc/c profile in Figure 2 would probably have been less negative at r < 0.3 Rq; an overall downward shift 
in the whole profile (such as results from the overall shift in the inferred helioseismic sound speed profile for 
smaller Rq values) would then probably lead to some improvement in rms{<5c/c}. However, the effect would 
still be much smaller than some of the other effects discussed below. 

Solar age effects: Figure 2 demonstrates that the uncertainty in the solar age i© has only a very minor 
effect on the solar sound speed profile — note that the shift illustrated here, from t© = 4.6 Gyr to 4.5 Gyr, is 
much larger than the observationally allowed range of solar ages, i.e., 4.55 Gyr < t© < 4.59 Gyr, as discussed 
in § 2. (These ages are defined to include the pre- main-sequence; main sequence ages can be obtained by 
subtracting 0.04 Gyr). The maximum allowed shift of 0.02 Gyr relative to the "best" solar age of 4.57 Gyr 
would yield negligibly small effects, namely, rms{Ac/c} w 0.0001 and rms{Ap/p} w 0.001, with maximum 
changes less than twice these values. Our results agree both qualitatively and quantitatively with the age 
sensitivity found in the recent work of Morel et al. (1997). 

Low-temperature opacity effects: Uncertainties in the low-temperature molecular opacities would not 
be expected to have much effect — in a convective region, such as the solar convective envelope, the structure 
is almost independent of the local opacity. As expected, using the Sharp (1992) molecular opacities ("«;sharp") 
below 10 4 K rather than the Alexander & Ferguson (1994) molecular opacities ( "^Alexander" ) led to essentially 
identical sound speed and density profiles — the "Ksharp" case is thus not plotted in Figure 3. The rms 
differences are negligible, less than a part in 10 4 for the sound speed and less than a part in 10 3 for the 
density. Only pre-main-sequence lithium depletion was significantly affected (see § 3.4). 
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Fig. 3. — Effect on the sound speed of opacity interpolation choices. The thin solid and dot-dashed 
lines show the effects of opacities where the C, N, O, and Ne abundances were respectively decreased 
("kopal:GN931C-No") or increased ("kopal:GN93TC-Nc") by their quoted errors of 15%, but the value of 
Z/X was not re-adjusted to reflect these changes, so these cases are not self-consistent (compare to the 
self-consistent cases in Fig. 5 below). Using opacity tables at a single constant metallicity ("const-Z K :" 
long-dashed line) leads to large opacity errors and thus large sound speed errors, which cannot be fixed by 
interpreting the metallicity error as "excess-CO" (short-dashed line). However, ignoring the effect on the 
opacity of relative differences in the abundance profiles of the individual metals results in only minor errors 
("approx-K:" dotted line). 

Interior opacity effects: Recently, Neuforge-Vcrhcccke ct al. (2001b) compared models using the 1995 
OPAL opacities (Iglesias & Rogers 1996; Rogers et al. 1996) with models using an updated version of the 
LEDCOP opacities from Los Alamos (Magee et al. 1995); under solar conditions, these two sets of opacities 
differ by up to 6% (just below the base of the convective envelope), although the authors indicate that about 
half of this difference is due to interpolation errors (from different temperature grids on which the opacities 
are tabulated). They find fractional sound speed differences up to Ac/c ~ 0.003 between solar models 
using these different opacity tables. Morel et al. (1997) have also demonstrated the serious impact of opacity 
changes on the sound speed and density profiles of solar models; they compared the 1995 OPAL opacities with 
the less- precise 1992 OPAL opacities (Rogers & Iglesias 1992) (albeit with models that neglected diffusion), 
finding that the improved opacities made an improvement of up to 0.005 in the sound speed and up to 0.03 
in the density (see their models SI and S2). Basu et al. (2000) compared a model with the 1995 OPAL 
opacities and the OPAL equation of state (Rogers et al. 1996) to a model with the 1992 OPAL opacities 
and the cruder Yale equation of state (Guenther et al. 1992) with the Debye-Hiickel correction (Bahcall, 
Bahcall, & Shaviv 1968) (their models did include diffusion); they likewise found an effect of up to 0.005 in 
the sound speed and up to 0.03 in the density, from the combination of these two changes in the input. We 
found that an even larger improvement of up to 0.007 in the sound speed resulted from changing from the 
even older Los Alamos (LAOL) opacities (Keady 1985, private communication) to the 1995 OPAL opacities, 
with an rms improvement of 0.004 as shown in Table 1 (there is also an improvement of up to 0.04 in the 
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density, with an rms improvement of 0.02) — however, such a large opacity change as this overestimates the 
uncertainty in the 1995 OPAL opacities (Rogers & Iglesias 1998). 

As pointed out by Morel et al. (1997), neglecting the opacity changes that result from metallicity 
variations in the Sun would lead to significant errors — e.g., errors of up to 0.0015 in the sound speed. In 
our reference standard solar model, we did our best to account for these temporal and spatial variations 
in the opacity due to these composition changes from diffusion and nuclear burning. As discussed in § 2, 
the metallicity value Z K that we used for metallicity interpolation in the OPAL opacity tables was scaled 
according to the changes in the elements heavier than oxygen ("Z K = Zh," where Zh is proportional to the 
heavy element abundance) . Several different methods were tested to account for the fact that changes in the 
CNO-element abundances (particularly in the solar core) are far from being proportional to changes in the 
heavy elements. As discussed in § 2, all such methods that we tested gave essentially identical results — even 
omitting the CNO-correction entirely had almost no effect, and only a negligible improvement resulted from 
a full "CNO-interpolation" case that included the full opacity effects of CNO abundance variations quite 
accurately (by interpolating among several separately-computed OPAL opacity tables with different CNO 
abundances). For our reference standard solar model, we simply assigned the sum of the non-proportional 
changes in the CNO abundances (namely, Z — Zh) to the "excess carbon and oxygen" interpolation variable 
of the OPAL opacity tables (i.e., GO ex = Z - Z h = Z - Z K ). 

An alternative approximation ("approx-K") is to set Z K ~Z (i.e., to interpolate the OPAL opacity tables 
in the local metallicity Z, but ignore effects of variations in the makeup of Z). As shown in Figure 3, this 
approximation yields results almost identical to those of our reference standard solar model; rms differences 
arc rms{Ac/c} « 0.0002 and rms{Ap/p} « 0.001, with maximum differences not very much larger. Morel 
ct al. (1997) compared two different ways of estimating the value of Z K (both being similar but not identical 
to our "approx-«" case); they likewise found only very minor differences in the sound speed between their 
two methods, but much larger effects on the density (see their models D3 and D12). 

We also tested cases where OPAL opacities had been calculated for mixes in which the abundances of 
C, N, O, and Ne were either all increased by their quoted errors of 15% ("kopal:GN93TC-Nc") or decreased 
by this amount ("kopal : gn931C-Nc"); relative to their abundances quoted by Grevesse & Noels (1993). 
Figure 3 illustrates an effect of up to 0.001 in the sound speed from such a change. However, we used a value 
of Z/X = 0.0245 for all of the models in Figure 3, which is not strictly consistent with such large abundance 
changes: since C, N, O, and Ne comprise the major portion of Z, a 15% change in their abundances should 
correspond to a change of ~ 12% in Z/X as well. We discuss such a self-consistent opacity-plus-composition- 
plus-Z/X variations under Z/X effects below. 

Even when using the most up-to-date OPAL opacities, one can still get significant errors if one neglects 
the effect on opacity of Z-changes (primarily due to diffusion) . The simplest case is to set Z K = Z ( "const- 
Z K "), where Zq is the protosolar metallicity; in effect, such a case uses only the OPAL opacity tables 
relevant to the protosolar metallicity and ignores the effect on the opacity of any subsequent changes in the 
metallicity. Figure 3 demonstrates that this "const- Z K " case yields errors of up to Ac/c ~ 0.0015 relative to 
the more accurate opacity interpolation of the reference standard solar model, in agreement with the results 
of Morel et al. (1997) (compare their models D10 and D12); the rms errors were rms{Ac/c} « 0.0010 and 
rms{Ap/p} « 0.007. One might attempt to fix up this neglect of metallicity variation by interpolation using 
the mildly CO-enhanced OPAL opacity tables, i.e., retaining a constant Z R = Z for opacity interpolation 
purposes but setting CO ex to the difference between Z K and the true value of Z that results from diffusion 
and nuclear burning, i.e., CO ea; = Z — Z K = Z — Z - Figure 3 illustrates that this approximation is 
a slight improvement over the "const-Z K " case but still not very satisfactory: it still has rms errors of 
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Fig. 4. — Effect on the sound speed of opacity uncertainties. The dot-dashed lines ("Iow-k" and "high-K") 
show the effect of a ~ 10% overall change in the opacities relative to the reference standard solar model 
(thick solid line), while the dashed lines show the effect of a ~ 5% overall change in the opacities. 



rms{Ac/c} « 0.0006 and rms{Ap/p} w 0.004. 

Interpolation errors can also arise from the finite grid spacing of the opacity tables in X, Z, T, and p. In 
creating the opacity interpolation routines, we tested the X- and Z-interpolation, finding that these should 
result in only minor errors (a fraction of a percent) in the opacity. The work of Neuforge-Verheecke et al. 
(2001b) suggests that T- and p-interpolation errors in the opacity can be larger, as much as a few percent. 

In addition to errors introduced by methods of interpolating in opacity tables, one must consider the 
errors in the actual opacity values contained in the tables. Such opacity errors will in general be functions of 
temperature and density; also, different elements will have different errors. Rogers & Iglesias (1998) estimate 
that there is a 4% uncertainty in the 1995 OPAL opacities from effects neglected in their calculations - 
Neuforge-Verheecke et al. (2001b) point out that the ~ 3% intrinsic differences between the OPAL and 
LEDCOP opacities are slightly less than this. Turcotte et al. (1998) showed that differences in the Rosseland 
mean opacities between 1992 and 1995 OPAL opacities do not exceed 7% over the run of temperature and 
density in the Sun's interior, and that these opacity differences yielded sound speed differences of up to 
0.002 in their solar models. In addition, errors in the observed relative heavy element abundances in the 
solar envelope will translate into opacity errors, since different elements have somewhat different opacities. 
Rogers & Iglesias (1998) estimate that such abundance uncertainties correspond to opacity uncertainties of 
order 5% at temperatures where ionization effects of the relevant elements yield a large contribution to the 
opacity (e.g., near 2 x 10 6 K for oxygen or neon). Turcotte et al. (1998) found that taking into account the 
opacity effects due to changes in the relative abundances of all the individual elements in Z led to opacities 
that differed by up to 2% from the opacities tabulated for the standard scaled-solar mctallicity, yielding 
sound speed differences of up to 0.001 and density differences up to 0.005 — these are several times larger 
than the effects discussed above that we found when testing effects of variations of C, N, and O relative to 
the heavier elements. 
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Fig. 5. — Effect on the sound speed of uncertainties in the observed solar surface composition. Relative 
to the reference standard solar model at Z/X = 0.0245 {thick solid line), switching to the more recent 
Z/X = 0.023 value of Grevesse & Sauval (1998) without including the corresponding changes in relative 
metal abundances (thin dotted line) has a larger effect than the full case with self-consistent abundances and 
opacities ("kopal:GS98 : " thin solid line)] the same is true when considering the older Z/X = 0.0277 cases of 
Grevesse (1984) (thin short-dashed line vs. self-consistent "KoPAL:Gr84" thick long-dashed line). Nonetheless, 
a self-consistent test of the uncertainties in the Grevesse & Sauval (1998) C, N, O, and Ne abundances (thick 
dot-dashed lines) shows a large effect (compare to the thin solid line). 



We did not attempt detailed element-by-element variations of the OPAL opacities in our models; nor 
did we test the effect of opacity variations in limited density or temperature ranges. Instead, we obtained 
a rough estimate of the maximum possible effects of uncertainties in heavy-element opacities by making 
an overall shift in the metallicity value Z K used for interpolation in the OPAL opacity tables. Figure 4 
illustrates the cases Z K = 0.9 Z ("low-/t"), Z K = 0.95 Z,Z K = Z ("approx-/t"), Z K = 1.05 Z, and Z K = 1.1 Z 
("high-/t"). The "low-/t" and "high-ft" cases correspond to an average shift in the opacities of order 10% over 
the solar interior relative to the "approx-ft" case (2 — 5% for r < 0.4 R Q , ~ 10% for 0.4 R Q < r < 0.7 Rq, 
~ 15% for 0.7 Rq < r < 0.92 Rq, and ~ 5% for r > 0.92 Rq). As may be seen from Figure 4, such an 
opacity shift of order 10% yields sound speed changes of up to Ac/c ~ 0.003, with rms{Ac/c} w 0.0016 
and rms{Ap/p} w 0.014. However, such a large opacity shift almost certainly overestimates the effects of 
opacity uncertainties. Rogers & Iglesias (1998) estimate that there is a 4% uncertainty in the 1995 OPAL 
opacities from effects neglected in their calculations, and it would be surprising if these yielded a uniform 
shift in the opacity throughout the Sun. Thus a better estimate of the effects of opacity uncertainties would 
be Ac/c - 0.001 and rms{Ap/p} - 0.005. 

Solar abundance (Z/X) effects: Our reference standard solar model used the observational value of 
Z/X = 0.0245 from Grevesse & Noels (1993), since their mixture was the one for which the standard OPAL 
opacity tables ("kopal:GN93") were available. The close-dashed curve in Figure 5 demonstrates the effects of 
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using the Grevesse & Noels (1993) relative metal abundances with the corresponding "kopal:GN93" OPAL 
opacities, but using a 13% higher value of Z/X — 0.0277, the older value that had been recommended 
by Grevesse (1984); the maximum sound speed difference relative to the reference standard solar model is 
Ac/c ~ 0.0030, with rms{Ac/c} « 0.0018 and rms{Ap/p} w 0.017. The dotted curve in Figure 5 illustrates 
a similar case with a 6% lower value of Z/X = 0.023, as recommended by the more recent work of Grevesse 
& Sauval (1998); it has maximum Ac/c - 0.0016, with rms{Ac/c} « 0.0010 and rms{Ap/p} « 0.009. (Note 
that most of the above effect comes from the different opacity that results from the changed solar Z value, 
as may be seen by comparing with the "high-K" and "1ow-k" curves in Fig. 3.) However, these comparisons 
are not strictly self-consistent, since it is the changes in the individual elemental abundances of the metals 
that add up to yield the changed Z/X ratio. Using the old abundance pattern of Grevesse (1984) and newly- 
computed OPAL opacities appropriate to it ("KoPAL:Gr84") leads to the wide-dashed Z/X = 0.0277 curve 
in Figure 5, with a slightly smaller maximum sound speed difference (of 0.0025), and rms{Ac/c} w 0.0014 
and rms{Ap/p} « 0.011. Similarly, using the Grevesse & Sauval (1998) abundance pattern and appropriate 
OPAL opacities ("kopal:GS9s") leads to the thin solid Z/X = 0.023 curve in Figure 5, reducing the maximum 
sound speed difference to 0.0012, with similarly reduced rms{Ac/c} « 0.0006 and rms{Ap/p} w 0.004. 

Neuforge-Verheecke et al. (2001a) performed a comparison identical to this last case, finding essentially 
the same effect on the sound speed, both qualitatively and quantitatively (maximum effect ~ 0.0018). Morel 
ct al. (1997) considered the effect of a 6% increase in Z/X, finding a maximum difference of 0.0007 in their 
sound speed and 0.003 in their density; this would imply a significantly lower sensitivity to Z/X than we 
found. This is probably due to the fact that the models in which they tested Z/X variations did not consider 
the effect on the opacities of the temporal and spatial variations in the heavy element abundances that arise 
from diffusion, but merely used opacities appropriate to a constant metallicity equal to the protosolar value 
(Z K = Z , as in our "const- Z K " case discussed above). 

Strictly, the uncertainty resulting from observational solar abundance errors can be estimated by varying 
the solar abundance values of Grevesse & Sauval (1998) within their quoted uncertainties, obtaining OPAL 
opacities with these revised compositions, calculating the resulting Z/X values, and running solar models 
with these self-consistent sets of input values. We have done this for two cases. Rather than performing large 
numbers of random variations of the abundances, we tested a case which should give something close to the 
maximum effect. The elements C, N, O, and Ne not only comprise the major part of the metallicity but also 
have relatively large errors of ~ 15%, and unlike other elements with large errors one cannot get a "better" 
value by using the meteoritic abundance instead. We therefore considered cases where C, N, O, and Ne were 
either all increased by 15% (Z/X = 0.0257, "k pal : gs98TC-Nc") or all decreased by 15% {Z/X = 0.0203, 
"ftOPAL:GS98iC-Ne") — i- e -i these self-consistent abundance variations correspond to 12% variations in Z/X. 
As illustrated in Figure 5, these cases lead to variations in the solar sound speed of up to 0.003 relative to 
the Z/X = 0.023 "k pal : gS98" case, with rms{Ac/c} « 0.0017 and rms{Ap/p} « 0.012. Of all the "input" 
uncertainties that we considered, these uncertainties in the solar abundances have the largest impact. 

Nuclear rate effects: Figure 6 demonstrates that the uncertainty in the 1 H(p, z/e+) 2 H nuclear burning 
rate (the pp rate) has a significant impact on the solar sound speed. Our reference standard solar model 
used the recommended nuclear rates from the NACRE compilation (Angulo ct al. 1999). These authors 
also supply "high" and "low" cases to indicate the allowed uncertainty range of each nuclear rate; in the 
case of the pp reaction, the high case is 8% above the recommended rate and the low case is 3% below 
it. We have tested the effects of nuclear rate uncertainties by computing variant standard solar models 
using high and low NACRE rate values. Figure 6 demonstrates that a high pp rate is preferable, if all 
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Fig. 6. — Effect on the sound speed of uncertainties in nuclear rates. The ~ 5% uncertainty in the basic pp 
rate (dot-dashed lines) has a much larger effect than the 20% uncertainty in the 3 He + a rate (dotted lines); 
other rate uncertainties have negligible effects, and are not plotted. The difference between the reference 
standard solar model using the NACRE rates (solid line) and the model using the Caughlan & Fowler (f 988) 
rates ("CF88:" dashed line) is largely due to the higher pp rate adopted by the latter authors. 



other parameters are kept constant: the high pp rate gives good agreement with the helioseismic reference 
profiles, except in the Sun's central regions where the helioseismic observations are the poorest. Our models 
indicate that a change of 5% in the pp rate yields changes of up to 0.003 in the sound speed, (0.0014 in the 
regions accurately probed by helioseismology, outside the core); the rms changes in such a case would be 
rms{Ac/c} « 0.0009 and rmsjAp/p} w 0.0f8. Antia & Chitrc (1999) also tested the effects of changes in 
the pp rate on the helioseismic profiles, concluding that a relatively high pp rate is preferred, consistent with 
our results discussed above. 

Figure 6 also demonstrates that the uncertainty of ±20% in the 3 He(a,7) 7 Be reaction leads to only a 
minor effect: a maximum sound speed change of 0.001 (or 0.0003 outside the core), with rms{Ac/c} w 0.0002 
and rms{Ap/p} « 0.004. Basu et al. (2000) considered the effect of setting the 3 He(a,7) 7 Be rate to zero; 
they found large effects from such an unphysically extreme change. Setting the rate to zero is equivalent 
to a 100% change, 5 times as large as the 20% change that we considered; thus it is consistent that their 
published effect is about 5 times as large as ours. 

We also tested the effects of the ±6% uncertainty in the 3 He( 3 He, 2p) 4 He reaction and of the ±30% 
uncertainty in the 14 N(p,j) 15 reaction (which determines the CNO-cycle rate). Such changes in these 
rates led to negligible effects on the sound speed and density profiles; we have not plotted these profiles in 
Figure 6, since they would be essentially superimposed on that of the reference standard solar model. 

ft is not surprising that the uncertainty in the pp rate has the largest effect on the sound speed and 
density profiles, since it is the basic rate that determines the overall p-p chain burning rate. 

We also computed a model using the previous standard set of nuclear rates, namely, the Kellogg nuclear 
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rate compilation of Caughlan & Fowler (1988). The resulting sound speed and density profiles are shown in 
Figure 6. With us still in Kellogg, carrying out this work in an office directly below his long-time office, it 
is especially gratifying for us to see that Willy Fowler's last published pp rate yields such good agreement 
with the current helioseismic reference profiles (the largest differences being near the center, where the 
observations are least accurate). 

Electron screening effects in nuclear rates: Gruzinov & Bahcall (1998) performed a careful quantum 
mechanical computation of the effects of electron screening on nuclear reactions in the Sun, and found that 
the Salpeter (1955) formula leads to only a very slight overestimate of the nuclear rates: 0.5% for the p + p 
reaction, 1.7% for the 3 He + 4 He reaction, 1.5% for the p + 7 Be reaction, and 0.8% for the p + 14 N reaction. 
Since these corrections are an order of magnitude smaller than the uncertainties given in the corresponding 
NACRE nuclear rates (Angulo et al. 1999), there was no point in considering separately the uncertainty in 
the screening corrections. Note that if one used the intermediate screening formulae (including the partial 
degeneracy correction) of Graboske et al. (1973), one would overestimate the nuclear reaction rates in the 
Sun by a further 1 to 3% relative to the weak screening formula — Gruzinov & Bahcall (1998) quoted a 
much larger effect in the opposite direction, but this was for a version of the Graboske et al. (1973) formulae 
that assumed completely degenerate electrons. 

Shaviv & Shaviv (1996) suggested that electron cloud-cloud interactions would increase the electron 
screening factor in the exponent of nuclear rates by a factor of 2/3; this was based on a "fundamental 
misconception concerning the dynamics of the interaction" (Briiggen & Gough 1997). Carraro, Schaier, & 
Koonin (1988) suggested that "dynamic screening" should reduce the Salpeter (1955) screening factor, since 
half of the screening effect comes from ions, which should not be able to adjust to the rapid motion of 
colliding/fusing nuclei. However, Gruzinov (1998) gave a general argument showing that in an equilibrium 
plasma there should be no such reduction, and Brown & Sawyer (1997) showed explicitly that such a 
modification of the Salpeter screening factor is exactly cancelled when one considers processes whereby 
Coulomb interactions with the colliding particles cause plasma excitations and de-excitations. Bahcall et 
al. (2002) point out that the more recent "dynamic screening" calculations of Shaviv & Shaviv (2000, 2001) 
seem likewise to be based on misconceptions. 

Bahcall et al. (2002) showed several ways of deriving the fact that the Salpeter (1955) electron screening 
formula gives the correct leading term under weak screening conditions, and used the rigorous formulation 
of Brown & Sawyer (1997) to show that, in the Sun, corrections to this leading term are small (of order 1%). 
They also pointed out that Tsytovich's alternative "anti-screening" formula (Tsytovich 2000; Tsytovich & 
Bornatici 2000), which yields a reduction in nuclear rates rather than an increase, yields unphysical results 
in two different limits. Nonetheless, some authors (Fiorentini, Ricci, & Villantc 2001; Weiss, Flaskamp, & 
Tsytovich 2001) have looked into the effects of using Tsytovich's alternative "anti-screening" formula, and 
have shown that it yields solar models that are not compatible with the helioseismic sound speed profile. 

Non-Maxwellian ion velocity distributions: Coraddu et al. (1999) suggest that the high-energy tail 
of the Maxwellian ion velocity distribution may be slightly depleted in the Sun; to first order, they express 
this modified distribution by f(E) ~ (fc7 1 )- 3 / 2 e -- E /' £ ' r - <5 (- E /' £T ) j where they estimate that 6 should be of 
the order 0.01. Such a large distortion of the high-energy tail would reduce the p + p reaction rate by more 
than 20% (and other pp-chain reaction rates by an order of magnitude); as may be seen by considering 
Figure 6, such a large reduction in the p + p reaction rate would not be compatible with the helioseismic 
sound speed profile. In addition, Bahcall et al. (2002) have criticized the above estimate of the magnitude 5 
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Fig. 7. Effect on the sound speed of uncertainties in diffusion constants. The tested cases of a 20% 
variation in the diffusion constant for helium (dot-dashed lines) may overestimate that uncertainty slightly; 
on the other hand, the tested cases of a 40% variation in the heavy element diffusion constants (dotted and 
short-dashed lines) may slightly underestimate the corresponding uncertainties — the reference standard 
solar model is shown by the solid line. 



of the effect, claiming that it should be negligibly small. Nonetheless, Turck-Chieze et al. (2001b) have 
considered the effect of such a Maxwellian distortion, with a smaller value of 5 = 0.002 — this reduces the 
p + p reaction rate by only 5%, yielding a solar model whose sound speed differs by 5c/ c < 0.005 from the 
hclioseismic profile (i.e., disfavored by helioseismology, but perhaps not ruled out entirely). 

Diffusion effects: There arc uncertainties in the diffusion coefficients; Proffitt (1994) estimates a 15% 
uncertainty in the diffusion constant of helium relative to hydrogen, and a <~ 50% uncertainty in the diffu- 
sion constant for oxygen relative to hydrogen. Figure 7 shows that an increase or decrease of 20% in the 
helium diffusion constants has only a modest effect: a maximum change of 0.001 in the sound speed, with 
rms{Ac/c} ~ 0.0008 and rms{Ap/p} « 0.007. The effect of increasing or decreasing the heavy element 
diffusion constants by 40% has a slightly smaller effect: a maximum change of 0.0006 in the sound speed, 
with rms{Ac/c} w 0.0004 and rms{Ap/ p} « 0.004. Note that including the effects of different diffusion 
rates for different heavy elements would have an effect less than half as large as this, as shown by the results 
of Turcotte & Christensen-Dalsgaard (1998) and Turcotte et al. (1998). 



3.2. Solar Convective Envelope Depth 

One of the key results of helioscismic observations is a highly precise value for the position R ce of the base 
of the solar surface convective region: Basu & Antia (1997) report a value of R ce — 0.713 ± 0.001 R@. Our 
reference standard solar model is in agreement with this value, having R ce = 0.7135 Rq (see Table 1), with 
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values negligibly different (R ce = 0.7133 or 0.7134 Rq) if we used the OPAL instead of the MHD equation 
of state at log p < —2. The new Grevesse & Sauval (1998) solar abundance observations (implying Z/X = 
0.023) yield a barely-consistent value of R ce = 0.7157 Rq; the maximum allowed abundance variations 
about this value (i.e., low and high Z/X values of 0.0203 and 0.0277, respectively) are inconsistent with the 
observed value (yielding R ce = 0.7209 and 0.7098 Rq, respectively). The solar age uncertainty of ±0.02 Gyr 
does not make a significant difference in R ce (see Table 1). As far as the uncertainties in the nuclear reaction 
rates are concerned, only the pp rate has a significant effect on R ce , of ±0.002 Rq. Uncertainties in molecular 
opacities do not have a significant effect on R ce , but changing from the old 1985 LAOL opacities to the 1995 
OPAL opacities does yield a large improvement (of 0.006 Rq); the remaining uncertainties in the 1995 OPAL 
opacities might thus be expected to have a small but possibly significant effect on R ce . The uncertainty in 
the diffusion constant for helium does have a significant influence (±0.003 Rq), but uncertainties in the 
diffusion constants for the heavy elements do not (effects < 0.001 Rq). 

Note that the cases favored by the sound speed profiles at the one- to two-sigma significance level (high 
opacities, high Z/X, or high pp rate) are disfavored by the observed R ce value at about the same significance 
level. 



3.3. Solar Helium Abundance 

Another key result of helioseismic observations is a fairly precise value for the present solar envelope 
helium mass fraction Y e (this value is lower than the Sun's initial helium abundance, due to diffusion). Infer- 
ring the solar helium abundance requires the use of a (theoretical) equation of state, as well as helioseismic 
frequency observations of modes that probe the solar convective region, particularly the He II ionization zone 
(Richard et al. 1998). 

Table 2 quotes helioseismic Y e values from eight recent papers — note that most authors calculate two 
separate values, obtained using the OPAL and MHD equations of state, respectively. Clearly, the systematic 
errors are much larger than the quoted internal errors from the helioseismic inversions, as shown by the 
relatively large differences between values obtained with these two different equations of state, and by the 
differences between determinations by different investigators. The OPAL equation of state is expected to 
be more accurate than the MHD equation of state over the bulk of the convective envelope, since the MHD 
equation of state was designed for use at log p < — 2 (which occurs in the Sun at r > 0.942 Rq); however, the 
He II ionization zone occurs further out (0.975 Rq < r < 0.985 Rq) in a region the MHD equation of state 
was specifically designed for, and where it may actually be more accurate than the OPAL equation of state 
(see, e.g., Richard et al. 1998; Basu et al. 1999; Di Mauro et al. 2002). Possible reasons for the relatively 
large scatter among Y e values of different investigators (even when using the same equation of state) are 
discussed by Di Mauro et al. (2002); they point out that observational frequencies for high-degree oscillation 
modes may suffer from significant systematic errors, and that the inversion formulae used to obtain Y e 
ignore several physical effects that may significantly affect the frequencies in this region of the Sun, such 
as nonadiabaticity, effects of mode excitations, or flows resulting from convective turbulence — certainly 
there is evidence that turbulent motions affect the /-modes (Gough 1993; Chitre, Christensen-Dalsgaard, & 
Thompson 1998; Murawski, Duvall, & Kosovichev 1998; Medrek, Murawski, & Roberts 1999; Di Mauro et 
al. 2002). 

One may perhaps ignore the extremely low value of Y e « 0.226 of Shibahashi, Hiremath, & Takata 
(1999), since it was obtained not from a direct inversion but by fitting to a sound speed profile (that had in 
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turn been obtained via a helioseismic inversion) ; relatively poor accuracy for this method is also suggested by 
the fact that the simultaneous determination by these authors of the depth of the solar convective envelope 
yielded R ce ss 0.718 Rq, quite far from helioseismic value of R ce = 0.713 ± 0.001 Rq (Basu & Antia 1997). 
The remaining Y e values from Table 2 have a mean of 0.245, a median value of 0.248, and a total range 
of 0.23 < Y e < 0.254, with most of the values lying in the range 0.24 < Y e < 0.25. These results can be 
conveniently summarized as an "observed" value of Y e = 0.245 ± 0.005 (where most of the uncertainty is due 
to systematic effects). 

Our theoretical reference standard solar model is in excellent agreement with this, having Y e — 0.2424 
independent of whether we used the OPAL or the MHD equation of state at log p < — 2. Low and high 
Z/X values (0.0203 and 0.0277, respectively) yield Y e values of 0.2396 and 0.2510, respectively — still quite 
acceptable. Similarly acceptable changes of ±0.005 in Y e result from a 4% uncertainty in the opacities, a 
15% uncertainty in helium diffusion constants, or a 50% uncertainty in heavy-element diffusion constants. 
Uncertainties in the solar age, luminosity, and radius and in the nuclear rates have only a negligible effect 
on Y e (see Table 1). 

3.4. Solar Lithium Abundance 

The present observed solar surface lithium abundance is log s( 7 Li) = 1.10 ± 0.10 as compared to the 
initial value of log s( 7 Li) = 3.31 ±0.04 obtained from meteorites (Grevesse & Sauval 1998), where log e( 7 Li) = 
\og(Nu/Nn) + 12 for number densities Nu and Ah of lithium and hydrogen, respectively The solar surface 
lithium depletion factor fu, relative to its initial value, is thus observed to be fu = 160 ± 40. Solar 
surface lithium can be depleted due to three causes: (1) lithium burning during the pre-main-sequence 
evolution, when the surface convection still reaches deeply into the interior; (2) rotationally induced mixing 
on the main sequence, which transports lithium down from the convective envelope to regions hot enough 
for lithium burning; (3) mass loss on the main sequence, which can cause the convective envelope to move 
inwards and engulf lithium-depleted regions. In this paper, we only consider the first of these, namely, 
the pre-main-sequence lithium destruction; rotational mixing is beyond the scope of this paper, and main 
sequence mass loss is discussed in the companion "Our Sun V" paper (Sackmann & Boothroyd 2002). Our 
reference standard solar model had a pre-main-sequence lithium depletion factor fu = 24, as shown in 
Table 1. This is a relatively large pre-main-sequence lithium depletion, leaving only a factor of ~ 7 to be 
accounted for by main sequence effects (2) and (3) above. Note that the recent models of Brun et al. (1999) 
also appear to find solar pre-main-sequence lithium depletion by a factor of order 10 (in good agreement 
with our results, considering the sensitivity of pre-main-sequence lithium depletion to the input physics). 

It is worth noting that, in the past, pre-main-sequence lithium depiction in the Sun has often been 
ignored or assumed to be negligible (which is not the case, as our present models show). Earlier models 
did in fact show relatively little pre-main-sequence solar lithium depletion, e.g., a depletion factor of 3 was 
reported by Sackmann et al. (1993). These low lithium depletion factors may have been caused partly by 
the use of older opacity tables in earlier models, but a major contributing factor was neglect of gravitational 
settling (diffusion) of helium and the heavy elements. These earlier models thus had a lower metallicity 
during the pre-main-sequence stage, matching the present solar surface metallicity, rather than being about 
10% higher (as diffusion models indicate). Because of the strong metallicity dependence of pre-main-sequence 
lithium depletion, such models yielded relatively small depletion factors. 

Pre-main-sequence lithium depletion depends quite sensitively on the structure of the solar models during 
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that stage of evolution. We give our results in terms of the solar lithium depletion factor fu, which is shown 
in Table 1. The zoning did not affect the lithium depletion factor significantly; nor did the uncertainties in 
the solar age, luminosity, and radius. Changes in the equation of state in the outermost regions can have 
a small effect (< 30%) on the lithium depletion factor: the "OPALeos-midT" and "OPALeos-lowT" cases 
had lithium depletion factors of fu = 17 and 19, respectively, as opposed to fu = 24 for the standard 
"OPALeos-hiT" case. 

The pre-main-sequence lithium depletion is extremely sensitive to both low-temperature and high- 
temperature opacities; use of the Sharp (1992) molecular opacities instead of the Alexander & Ferguson 
(1994) ones halved the lithium depletion factor (fu = 10), as did our "1ow-k" test case (fu = 10), while our 
"high-/t" test case nearly tripled it (fu = 71). There is also a relatively large sensitivity to the uncertainty in 
the observed solar abundances. A low Z/X ratio of 0.0203 (with "kopal:GS98|C-Nc") yields only two-thirds 
as much lithium depletion (fu = 15), while a high Z/X ratio of 0.0257 (with "kopal:GS98TC-No") yields 
half again as much lithium depletion (fu = 35). 

Uncertainties in the diffusion constants can affect the lithium depletion factor significantly, primarily 
due to the effect of the different initial composition; our diffusion test cases have depletion factors ranging 
from fu = f 8 to 33. 

Except for the 7 Li + p rate, uncertainties in the nuclear rates have almost no effect on the extent of 
lithium depletion. For the 7 Li + p rate, the ±14% uncertainty quoted by the NACRE compilation (Angulo 
ct al. 1999) corresponds to an uncertainty of about ±50% in the depletion factor (i.e., a range in the from 
fu = 16 to 38). 

Solar beryllium abundance: The observed solar beryllium abundance is log e( 9 Be) = 1.40±0.09, consis- 
tent with no depletion relative to the meteoritic value of log e( 9 Bc) = 1.42 ± 0.04. The quoted uncertainties 
of these values imply that solar beryllium cannot have been depleted by more than a factor of 2 (3-<r upper 
limit). Our solar models all have negligible amounts of beryllium depletion, of order 1% — in good agreement 
with these observational beryllium results. Other authors have shown that parameterized rotational- mixing 
models which match the observed solar lithium depiction yield relatively minor beryllium depletion, in agree- 
ment with observations. For example, the rotational-mixing solar models of Brun et al. (1999), with total 
lithium depletion of a factor of order 100 (a factor of ~ 10 occurring on the main sequence), deplete beryllium 
by only about 12%. Note that rotational-mixing models which ignore pre-main-sequence lithium depletion 
(i.e., which over-estimate the main-sequence lithium depletion by an order of magnitude) might be expected 
to overestimate the beryllium depletion as well: for example, Richard et al. (1996) required main-sequence 
lithium depletion by a factor of 155 in their solar models with rotational mixing, and found that these models 
then implied beryllium depletion by a factor of 2.9. 



3.5. Solar neutrinos 

We will not devote much space to the predicted solar neutrino values, since it has long been concluded 
that matching the observed neutrino capture rates requires not revised astrophysics but new neutrino physics, 
e.g., Mikheyev-Smirnov-Wolfenstein (MSW) neutrino oscillation effects — even non-standard solar models 
(e.g., with core mixing or a low-metallicity core) cannot simultaneously satisfy the neutrino observations and 
the helioseismic constraints (see, e.g., Bahcall, Basu, & Pinsonneault 1998a; Suzuki 1998; Basu et al. 2000; 
Bahcall et al. 2001; Watanabe & Shibahashi 2001; Turck-Chiezc et al. 2001a,b; Choubey et al. 2001; Guzik 
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ct al. 2001). In Table 1 we present the theoretically predicted neutrino capture rates for the 37 C1 and 71 Ga 
experiments, and the predicted flux of 8 B neutrinos. As is normally obtained, our theoretical predicted 
neutrino rates are much in excess of the observed values, i.e., 6.4 to 8.9 SNU is predicted for the 37 C1 
experiment, as compared to the observed value of 2.56 ± 0.23 SNU (Davis 1994; Cleveland et al. 1998), and 
127 to 141 SNU for the 71 Ga experiments, as compared to the observed value of 74.5 ± 5.7 SNU (combined 
value from SAGE and GALLEX+GNO: Hampel et al. 1999; Abdurashitov et al. 1999; Altmann ct al. 1999; 
Gavrin 2001; Ferrari 2001). Likewise, the models predict 8 B neutrino fluxes of 4.4 to 6.3 x 10 6 cm~ 2 s _1 , 
as compared to the value of (2.32 ± 0.08) x 10 6 cm" 2 s -1 implied by the Super-Kamiokande measurement 
(Fukuda et al. 2001). 

4. Conclusion 

Hclioseismic frequency observations enable the adiabatic sound speed c and adiabatic index Ti to be 
inferred with an accuracy of a few parts in 10 4 , and the density p with an accuracy of a few parts in 10 3 . 
These quantities can also be computed on purely theoretical grounds. It is important to understand the 
uncertainties in these theoretical quantities (arising from uncertainties in the input data), when comparing 
them to the values inferred from helioseismic measurements. These uncertainties in the theoretical standard 
solar model are presented below. 

(1) Abundances of the elements: For the standard solar model, we found that the largest impact on 
the sound speed arises from the observational uncertainties in the photosphcric abundances of the elements. 
The key elements for which accurate meteoritic determinations are not available are C, N, O, and Ne, with 
uncertainties of 15% (leading to an uncertainty of order 10% in the solar Z/X ratio). We determined that 
this abundance uncertainty affects the sound speed profile in the solar model at the level of 3 parts in 10 3 . 

(2) OPAL opacities, pp nuclear rate, and diffusion constants: The estimated 4% uncertainty 
in the OPAL opacities, the ~ 5% uncertainty in the basic pp nuclear reaction rate, the ~ 15% uncertainty 
in the diffusion constants for the gravitational settling of helium, and the <~ 50% uncertainties in diffusion 
constants for the heavier elements, all affect the sound speed at the level of 1 part in 10 3 . 

(3) Solar radius and low-temperature equation of state: Different observational methods yield 
values for the solar radius differing by as much as 7 parts in 10 4 ; this leads to uncertainties of a few parts 
in 10 3 in the sound speed in the solar convective envelope, but has negligible effect on the interior (recall, 
however, that while the sound speed in the solar model is not affected, there is a systematic effect on 
the "observed" helioseismic sound speed profile). Uncertainties in the low-temperature equation of state 
(log T < 5.5) lead to uncertainties of order a part in 10 3 in both the sound speed and the adiabatic index Ti 
in the convective envelope. 

(4) Rotational mixing and high-temperature equation of state: We did not explicitly consider 
the effects of rotational mixing or uncertainties in the interior equation of state, but other investigators have 
found these to yield uncertainties in the sound speed of order a part in 10 3 , as discussed in § 3.1. 

(5) Other sources of uncertainty: We found that other current uncertainties, namely, in the solar age 
and luminosity, in nuclear rates other than the pp reaction and in the low-temperature molecular opacities, 
have no significant effect on the quantities that can be inferred from helioseismic observations (although 
some of these can have significant effects on neutrino fluxes and/or the extent of pre-main-sequence lithium 
depletion). 
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(6) Depth of envelope convection: Our reference standard solar model (with Z/X = 0.0245) 
yielded a convective envelope position R ce = 0.7135 Rq, in excellent agreement with the observed value of 
0.713 ± 0.001 Rq, and was significantly affected (±0.003 Rq) only by uncertainties in Z/X, opacities, the 
pp rate, and helium diffusion constants. 

(7) Envelope helium abundance: Our reference model yielded an envelope helium abundance 
Y e = 0.2424, in good agreement with the range of values inferred from helioseismic observations (which 
we summarize as a helioseismic Y e = 0.245 ± 0.005); only extreme variations in Z/X, opacities, or diffusion 
constants yielded Y e variations as large as 0.005. 

(8) Pre-main-sequence lithium depletion: For the standard solar model, the predicted pre-main- 
sequence lithium depletion is a factor of order 20 (an order of magnitude larger than that predicted by earlier 
models that neglected gravitational settling and used older opacities) . The lithium depletion factor can vary 
between ~ 10 and ~ 40 when one varies the input physics, i.e., it is uncertain by a factor of 2. 

(9) Solar neutrinos: For the standard solar model, the predicted neutrino capture rate is uncertain 
by ~ 30% for the 37 C1 experiment and by ~ 3% for the 71 Ga experiments (not including uncertainties in the 
capture cross sections), while the 8 B neutrino flux is uncertain by ~ 30%. 
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Tabic 1. Characteristics of Our Solar Models 3, 
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4 


.00008 


.00115 


24.06 
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.00049 
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.00168 


24.30 
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1.813 


.02029 
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.2425 
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.00091 
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24.36 
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135.2 


5.89 
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.00091 
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4 
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24.36 


7.26 


132.5 


4.77 
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high p + Li rate 


1.814 


.02030 
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133.8 


5.33 



Tabic 1 — Continued 
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45. low p + 7 Li rate 


1.814 


.02030 
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.2424 


.7136 


.00139 


.00091 


.01772 


4 


.00001 


.00005 


15.69 


7.89 


133.8 
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46. CF88 nuclear rates 


1.832 


.02018 


.2767 


.2440 


.7117 


.00077 


.00035 


.00500 


4 


.00090 


.01709 


16.28 


7.85 


130.4 


5.46 


47. 20% low D Ho 


1.781 


.02018 


.2760 


.2485 


.7169 


.00217 


.00162 


.02439 


4 


.00084 


.00682 


20.15 


7.71 


132.9 


5.20 


48. 20% high D Ha 


1.846 


.02038 


.2760 


.2368 


.7106 


.00080 


.00046 


.01171 


4 


.00081 


.00635 


29.03 


8.12 


135.2 


5.51 


49. 40% low D z 


1.824 


.01949 


.2725 


.2385 


.7128 


.00165 


.00130 


.02187 


4 


.00038 


.00420 


18.43 


7.42 


131.2 


4.98 


50. 40% high D z 


1.804 


.02110 


.2795 


.2463 


.7145 


.00129 


.00059 


.01412 


4 


.00033 


.00369 


32.46 


8.42 


136.9 


5.72 



a Mixing length parameter a, prc-solar mctallicity Zq and helium mass fraction Yq, present envelope helium abundance Y e , position R ce of the base of envelope 
convection, rms fractional sound speed and density differences relative to the Sun's inferred helioscismic profiles and relative to the reference standard solar 
model, pre-main-sequence lithium depletion factor / Li , predicted capture rates (in SNU) <I>ci an d ^Ga f° r chlorine and gallium experiments, respectively, and 
predicted flux $b of 8 B neutrinos fin units of 10 6 cm s — 1 ). 

b Reference standard solar model: OPAL EOS at log p > —1.5, MHD EOS at log p < — 2, high-T opacities kopal:GN93 interpolated in Z K — = 
iYl heavy ] / E heav y (^"* )°] as wcn as m "excess" C and O (such that C ex + O ex = CO ex — Z — Z K ), low-T opacities ^Alexander; NACRE nuclear rates, 
gravitational settling of He and heavy elements, Z/X — 0.0245, Lq — 3.854 x 10 33 erg s -1 , Rq — 695.98 Mm, and £q — 4.6 Gyr; both fine-zoned and 
coarse-zoned cases were computed. Variant models have the same input values, except as specified in the first column. 

c OPALcos-midT: use OPAL EOS at log T > 4.0, MHD EOS at log T < 3.9. 

d OPALcos-lowT: use OPAL EOS at log T > 3.75, MHD EOS at log T < 3.7 (on the main sequence, MHD EOS is only used outside the photosphere). 

° This reference standard solar model and all subsequent models in the table use the coarse zoning. 

f For these cases, the relative rms values compare the coarse-zoned cases to the corresponding fine-zoned cases. 

g For these cases, the relative rms values compare the variant models with the case relative to which the parameter variation was performed. 
h Most recent solar luminosity value of Lq — L oes t = 3.842 x 10 33 erg s _1 . 

1 OPAL opacities not interpolated in il cxcess" C and O, i.e., opacities computed as if the CNO clement abundance profiles in the Sun were always in the same 
proportion to the Fe profile (Z^ — Z^ ^ Z, but no correction term: CO ex — 0.0). 

J "Unmatched re" means that Z/X has been varied relative to the recommended value of the given mixture (for which the OPAL opacities had been computed). 
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Table 2. Helioseismically Measured Present Solar Envelope Helium Abundance Y e 



Y e Using OPAL Eos Y e Using MHD Eos Reference 





0.242 ± 0.003 


1 


0.249 ±0.001 


« 0.25 


2 


w 0.23 a 


« 0.25 a 


3 


0.248 ±0.001 


0.232 ±0.006 


4 


0.248 ±0.001 


w 0.252 


5 


0.248 ±0.002 


w 0.242 


6 


w 0.226 b 




7 


0.2539 ±0.0005 


0.2457 ±0.0005 


8 



a Obtained via an entropy calibration. 

b Obtained from a fit to a solar sound speed profile, not a 
direct helioseismic inversion. 
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